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ABSTRACT 

In  this  thesis,  a  code  division  multiple  access  (CDMA) 
communication  system  for  the  low  frequency  (LF)  channel  is 
proposed,  discussed  and  analyzed.  This  LF/CDMA  scheme  is  similar 
to  classical  CDMA  schemes  in  that  K  users  share  a  channel  by 
phase  modulating  their  transmissions  with  signature  sequences. 
Our  LF/CDMA  scheme  is  different  in  that  each  user's  signature 
sequence  set  consists  of  M  orthogonal  sequences  and  thus  log2M 
bits  of  information  are  transmitted  by  choosing  among  the 
signature  sequences.  Additionally,  the  users  use  r-phase 
modulation  and  our  model  includes  an  impulsive  ( non-Gauss i an ) 
noise  source  to  model  LF  atmospheric  noise. 

We  derive  a  locally  optimum  (small  signal)  receiver 
structure  for  our  LF/CDMA  scheme.  This  receiver  consists  of  a 
bandpass  correlator  followed  by  a  sampler,  a  zero  memory 
nonlinearity  and  M  discrete  time  matched  f ilter/cor relators.  We 
analyze  the  performance  of  this  structure  in  combined  multiple 
access,  impulsive  and  Gaussian  noise.  When  the  noise  is 
dominated  by  either  its  multiple  access  or  Gaussian  component, 
the  receiver  predominantly  operates  in  the  linear  region  of  the 


nonlinearity  and  performance  is  similar  to  that  of  a  linear 
receiver . 

We  analyze  the  multiuser  error  probability  (MEP)  of  the 
linear  receiver  by  using  a  Gaussian  assumption  to  find  an 
approximation  and  a  characteristic  function  method  to  achieve  a 
sharp  upper  bound.  The  approximation  and  the  bound  sometimes 
have  substantial  discrepancies  and  these  differences  are 
discussed. 

Finally,  we  design  two  actual  sequence  sets  for  our  LF/CDMA 
scheme.  Our  sequence  set  design  strategy  is  to  minimize  the 
maximum  magnitude  of  the  interference  that  any  one  user's 
transmitter  can  introduce  into  the  receiver  of  another  user.  Two 
sequence  designs  result,  where  one  is  based  on  additive 
characters  and  the  other  on  multiplicative  characters.  These 
designs  are  analyzed  and  their  differences  are  discussed. 
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CHAPTER  1 
INTRODUCTION 


1.1  Overview 

In  this  thesis,  we  design  and  analyze  a  code  division 
multiple  access  (CDMA)  communication  system,  which  is  especially 
suited  for  the  low  frequency  (LF)  band.  The  LF  band  is  the 
portion  of  the  frequency  spectrum  from  30  kHz  to  300  kHz.  General 
descriptions  of  propagation  modes,  noise  models  and  engineering 
practice  applicable  to  this  band  are  available  in  Naval  Shore 
Electronics  (1972)  and  Watt  (1967).  The  LF  band  is  characterized 
by  propagation  modes,  which  make  long  range  communication 
possible  and  by  impulsive  atmospheric  noise,  which  is  distinctly 
non-Gaussian  in  nature.  Most  importantly,  LF  communication 
channels  have  inherently  narrow  bandwidth  and  this  limitation  is 
a  most  important  consideration  when  designing  LF  spread  spectrum 
systems.  The  system,  which  we  discuss  and  analyze  in  this 
thesis,  is  designed  to  make  thorough  use  of  this  limited 
bandwidth.  As  such,  it  has  many  differences  (as  well  as  many 
similarities)  with  the  classical  spread-spectrum  communication 
schemes,  such  as  those  described  in  Pursley  (1977,1981);  Pursley 
and  Sarwate  (1977);  Pursley,  Garber  and  Lehnert  (1980);  and 
Garber  and  Pursley  (1981). 

Our  low  frequency  code  division  multiple  access  (LF/CDMA) 
communication  system  is  similar  to  the  classical  model  in  the 
following  respects.  First  of  all,  it  has  K  users  sharing  a 
channel  by  phase  modulating  their  transmissions  with  signature 
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sequences.  These  users  make  no  attempt  at  frequency  separation 
and  the  scheme  is  asynchronous,  which  means  that  the  users  are 
not  time  coordinated  in  any  way.  The  model  also  includes  an 
additive  noise  source  to  account  for  atmospheric  and  receiver 
noise . 

Our  scheme  also  has  the  following  features,  « a i c h  are  not 
part  of  the  classical  models.  Each  user's  signatu  .  iquence  set 
consists  of  M  orthogonal  sequences  and  thus  log2M  bits  of 
information  are  transmitted  by  choosing  among  the  signature 
sequences.  Typically,  M  is  of  the  order  of  several  hundreds. 
Each  signature  sequence  is  r-valued  (which  means  r-phase 
modulation  is  used),  where  r  may  be  as  much  as  one  hundred.  We 
can  consider  size  M  sequence  sets  and  r-phase  modulation,  because 
the  long  times  associated  with  the  LF  band  allow  an  increased 
freedom  in  signal  set  design.  Furthermore,  the  initial  carrier 
phase  for  each  user  is  not  modelled  as  random,  because  LF 
transmitter  techniques  are  such  that  the  carrier  phase  is  fixed 
and  stable  with  respect  to  the  signal  envelope.  The  net  result 
is  that  each  user  produces  transmitted  signals  with  bandwidths  of 
5  kHz  to  20  kHz  and  with  information  rates  up  to  1  kilobit  per 
second . 

As  mentioned  above,  LF  atmospheric  noise  is  "impulsive",  and 
not  Gaussian.  We  discuss  the  various  models,  which  have  been 
used  to  describe  the  first  order  probability  function  (pdf)  of 
impulsive  noise  and  present  an  argument  for  using  a  Cauchy  pdf. 
After  discussing  atmospheric  noise  models,  we  derive  a  locally 


optimum  Bayes  detector  for  use  as  a  receiver.  This  receiver 
contains  a  filter  matched  to  the  basic  "chip"  waveform  followed 
by  a  zero  memory  nonlinearity  (ZNL).  The  ZNL  is  in  turn  followed 
by  a  bank  of  M  correlators,  where  each  correlator  corresponds  to 
one  of  the  possible  M  signature  sequences.  We  analyze  t  ne 
performance  of  our  receiver  and  find  that  the  performance  is  a 
strong  function  of  Fisher's  information  number,  which  completely 
depends  on  the  statistics  of  the  filtered  noise.  The  analysis 
also  shows  that  all  LF/CDMA  receivers  should  employ  a  ZNL  to 
limit  the  effect  of  the  impulsive  noise.  Additionally,  if  the 
multiple  access  or  Gaussian  noise  is  strong,  then  the  performance 
of  this  nonlinear  receiver  closely  resembles  the  performance  of  a 
linear  receiver  operating  in  the  combined  multiple  access  and 
Gaussian  noise  (without  impulsive  noise.) 

Consequently,  we  consider  a  linear  receiver  (without 
impulsive  noise),  which  we  analyze  using  a  characteristic 
function  method  (Geraniotis  and  Pursley  (1982)).  Tne 
characteristic  function  approach  yields  a  bound  for  tne 
receiver's  probability  of  error.  An  error  probability 
approx imat  ion  can  be  obtained  by  assuming  the  multiple-access 
interference  is  approximately  Gaussian  and  using  the  well  known 
results  for  coherent  M-ary  communication  in  Gaussian  noise.  The 
results  from  the  characteristic  function  analysis  can  differ 
quite  significantly  from  those  based  on  the  Gaussian  assumption, 
and  we  provide  a  brief  discussion  of  these  discrepancies.  For 
most  of  our  linear  analysis,  we  consider  the  user  codes  to  o> 
random  cosets  of  an  orthogonal  code. 


Following  the  system  analysis  portion  of  the  thesis,  we 
design  two  specific  sequence  sets,  which  are  well  suited  for  our 
LF/CDMA  communication  scheme.  Both  designs  have  upper  bounds  for 
the  maximum  value  of  the  magnitude  of  the  multiple  access 
interference.  These  interference  bounds  are  determined  through 
the  use  of  a  partial  sum  theorem.  Even  though  our  main  interest 
in  this  thesis  is  multiple  access  interference,  skywave 
(specular  multipath)  and  signal  acquisition  performance  could  be 
major  issues  in  the  design  of  any  LF  communication  systems. 
Consequently,  both  of  our  sequence  sets  are  also  designed  to  have 
good  skywave  and  acquisition  performance. 

One  of  our  sequence  set  designs  is  based  on  additive 
characters  and  the  other  is  based  on  multiplicative  characters. 
The  main  differences  between  the  two  designs  is  the  number  of 
required  phases  and  number  of  accommodated  users.  The  additive 
character  design  accommodates  N  users,  but  requires  N  phases, 
where  N  is  the  sequence  length.  The  multiplicative  character 
design  uses  many  fewer  than  N  phases,  but  consequently 
accommodates  many  fewer  than  N  users. 

In  the  remainder  of  this  introductory  chapter,  we  will 
describe  spread  spectrum  multiple  access  research,  the  thesis 
outline,  and  the  contributions  of  the  thesis. 


1.2  Spread  Spectrum  Multiple  Access  Communication 


Spread  spectrum  systems  are  those  where  the  transmitted 
signal  occupies  much  greater  bandwidth  than  the  data  signal 


itself.  Bandspreading  is  commonly  achieved  either  by  frequency 
hopping  (FH),  where  the  center  frequency  of  the  transmitted 
signal  is  varied,  or  direct  sequence  (DS)  which  varies  the  phase. 
Spread  spectrum  systems  use  greater  bandwidth,  because  the 
transmitted  signal  has  been  carefully  designed  to  achieve 
specialized  system  goals,  some  of  which  are: 

•  high  resolution  ranging 

•  multipath  rejection 

•  signal  hiding 

•  multiple  access  capability. 

The  last  of  these  goals  is  of  greatest  interest  to  us  and  we  now 
consider  a  direct  sequence  spread  spectrum  multiple  access 
system. 

Some  spread  spectrum  systems  use  phase  or  direct  sequence 
modulation  to  achieve  multiple  access  communication.  Such 
systems  are  called  code  division  multiple  access  (CDMA)  systems 
and  one  such  scheme  has  been  modelled  and  analyzed  by  Pursley 
(1977).  In  Pursley's  model,  each  transmitter  has  a  single  binary 
signature  sequence  and  phase  modulates  its  carrier  with  the 
sequence  or  its  negative  depending  on  the  data.  Each  receiver 
contains  a  filter  matched  to  the  carrier  modulated  by  the  desired 
sequence  and  uses  the  output  of  the  filter  as  its  decision 
statistic.  The  receiver  decides  the  value  of  the  data  bit  based 
on  the  sign  of  the  decision  statistic.  However,  the  receiver's 
decision  is  not  perfectly  reliable,  because  of  noise  in  the 
channel  (which  is  assumed  to  be  an  independent  additive  white 
Gaussian  process)  and  because  of  signals  from  each  of  the 


competing  users,  who  share  the  channel.  Each  of  these  users 
contributes  an  interference  term  to  the  decision  statistic  and 
the  sum  of  these  terms  is  known  as  multiple  access  interference. 
Each  of  the  interference  terms  is  a  function  of  the  following 
var iables : 


•  the  time  of  arrival  delay  of  the  competing  users' 
signal  with  respect  to  the  desired  signal 

•  the  phase  offset  of  the  competing  users' 
transmitters  relative  to  the  desired  signals' 
transmitter 

•  the  competing  users'  data. 

In  general,  these  variables  may  be  characterized  in  dif¬ 
ferent  ways  to  fit  the  system  under  consideration  and  the  issues 
at  hand.  In  Pursley's  1977  paper,  two  different  analyses  based 
on  different  characterizations  of  the  underlying  variables  are 
presented  and  are  quickly  reviewed  below. 

Worst  case  analysis  considers  the  delay,  phase  offset,  and 
data  variables  to  be  deterministic  and  finds  the  value  of  these 
variables  for  which  the  multiple  access  interference  takes  its 
maximum  magnitude.  This  approach  results  in  an  upper  bound  on 
the  receiver's  probability  of  error  and  provides  a  mini-max 
design  criterion  for  the  selection  of  signature  sequences.  This 
mini-max  sequence  design  criterion  is  further  investigated  in 
Pursley  and  Sarwate  (1977),  Sarwate  (1979),  and  Sarwate  and 
Pursley  (1980).  The  deterministic  analysis  also  shows  that  the 
worst  case  performance  is  independent  of  the  shape  of  the  chip 


waveform  (the  chip  waveform  is  the  basic  signal,  which  the 
transmitter  phase  modulates  and  transmits).  Worst  case  analyses 
guarantee  a  performance  minimum  under  the  worst  of  conditions. 
On  the  other  hand,  the  worst  of  conditions  generally  occurs  so 
infrequently  that  no  estimate  of  the  average  performance  of  the 
system  is  obtained. 

For  this  reason,  average  system  performance  is  estimated  by 
considering  each  interference  term  to  be  a  function  of 
independent  random  variables.  This  approach  implies  that  the 
multiple  access  interference  is  also  a  random  variable  and  the 
mean  and  variance  of  the  multiple  access  interference  may  be 
calculated.  The  mean  is  zero  and  the  variance  may  be  used  to 
calculate  the  receiver's  signal  to  noise  ratio.  The  signal  to 
noise  ratio  may,  in  turn,  be  used  to  calculate  the  receiver's 
probability  of  error  if  the  multiple  access  interference  is 
assumed  to  have  a  Gaussian  distribution.  This  Gaussian 
assumption  gives  reasonably  accurate  estimates  of  the  actual 
probability  of  error  (Pursley  et  al.  (1982),  Geraniotis  and 
Pursley  (1982)).  In  Pursley's  model,  the  average  performance 
analysis  yields  signal  to  noise  ratios  (and  consequently 
probability  of  error  estimates)  as  a  function  of  the  mean  square 
correlation  parameters  of  the  signal  sets.  These  parameters  are 
explored  as  alternate  sequence  design  criteria  in  Pursley  and 
Sarwate  (1977)  and  Sarwate  and  Pursley  (1980).  Additionally,  the 
average  analysis  shows  that  the  signal  to  noise  ratio  is  a 
function  of  the  mean  square  correlation  parameters  of  the  chip 
waveforms.  These  parameters  have  been  used  as  chip  waveform 


design  guides  by  Lehnert  (1981). 

Random  sequence  analysis  is  also  explored  as  a  CDMA  system 
design  tool  by  Pursley  (1977)  and  Roefs  and  Pursley  (1976). 
Here,  the  K  signature  sequences  are  considered  to  be  mutually 
independent  random  sequences  and  each  sequence  is  considered  to 
be  a  sequence  of  independent  identically  distributed  binary 
random  variables.  The  expected  value  of  the  mean  square 
correlation  parameters  for  the  random  sequences  may  then  be 
readily  calculated.  In  this  way,  random  sequence  analysis  pro¬ 
vides  a  straight  forward  estimate  of  system  performance  (signal 
to  noise  ratio  and  probability  of  error)  without  requiring  any 
actual  sequence  design. 

Recently,  the  average  probability  of  error  in  CDMA  systems 
has  been  calculated  via  a  characteristic  function  method 
(Geraniotis  and  Pursley  (1982)  and  Geraniotis  (1982)).  This 
approach  eliminates  the  uncertainty  due  to  the  Gaussian 
approximation  and  has  shown  that  the  latter  does  sometimes  lead 
to  appreciable  errors.  In  Geraniotis  and  Pursley  (1982),  the 
characteristic  function  method  is  applied  to  binary  and 
quaternary  CDMA  systems,  operating  in  additive  white  Gaussian 
noise.  Results  are  given  for  various  signature  sequences  (m- 
sequences  and  Gold  sequences)  and  a  collection  of  chip  waveforms. 
In  Geraniotis  (1982),  the  method  is  used  on  direct  sequence  and 
frequency  hopped  SSMA  communication  systems,  which  suffer  fading 
as  well  as  additive  noise.  In  this  latter  work,  specific 
signature  sequences  are  not  considered,  because  random  sequence 
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analysis  is  used. 

The  role  of  all  these  CDMA  design  tools  in  this  thesis  will 
be  discussed  shortly. 

1.3  Thesis  Outline 

In  Chapter  2  of  the  thesis,  we  describe  atmospheric  noise 
and  consider  appropriate  nonlinear  receivers.  Specifically  in 
Section  2.2,  we  introduce  our  LF/CDMA  model  and  in  Section  2.3, 
we  describe  LF  atmospheric  noise  and  our  choice  for  modelling  its 
first  order  probability  distribution  function.  In  Section  2.4, 
we  derive  a  locally  optimum  Bayes'  detector  ( LOBD )  of  M  signals 
and  in  Section  2.5  we  approximate  the  performance  of  this 
structure  in  combined  atmospheric,  Gaussian,  and  multiple  access 
noise . 

In  Chapter  3,  we  consider  linear  receivers  and  calculate 
average  multiple  access  performance.  In  Section  3.2,  we  begin  by 
considering  single  user  systems  and  in  Section  3.3,  we  present 
the  random  coding  ideas  which  we  will  use  in  our  multiuser 
analysis.  Section  3.4  is  devoted  to  multiuser  analysis  and  both 
the  Gaussian  approximation  and  the  characteristic  function  method 
are  used.  This  section  compares  the  results  of  these  two 
approaches  and  also  considers  some  alternate  signalling  schemes. 

In  Chapter  4,  we  design  specific  sequence  sets  for  our 
LF/CDMA  scheme  by  using  a  minimax  approach.  In  Section  4.2,  we 
begin  by  developing  a  minimax  sequence  design  criterion  by 
performing  a  worst  case  multiple  access  interference  analysis. 
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The  next  section  develops  this  criterion  through  application  of  a 
partial  sum  theorem.  In  Sections  4.4  and  4.5,  we  find  additional 
sequence  design  criteria  for  multipath  and  signal  acquisition 
respectively.  In  Section  4.6,  we  describe  additive  characters  in 
general  and  use  them  to  design  sequences  for  our  LF/CDMA  scheme. 
In  Section  4.7,  we  describe  multiplicative  characters  and  use 
them  to  design  LF/CDMA  sequences. 

In  Chapter  5,  we  summarize  the  main  conclusions  of  the 
thesis . 


1.4  Contributions  of  the  Thesis 

The  first  contribution  of  the  thesis  is  the  low  frequency 
CDMA  signalling  scheme  itself.  Specifically,  each  user  has  M 
orthogonal  signature  sequences  and  consequently  log2M  bits  of 
information  are  communicated  by  choosing  among  the  sequences. 
Additionally,  the  sequences  are  r-valued,  which  means  that  each 
user  is  employing  r-phase  modulation. 

The  second  major  contribution  of  the  thesis  is  the 
consideration  of  impulsive  atmospheric  noise.  This  noise  causes 
us  to  derive  a  locally  optimum  Bayes'  detector  of  M  signals, 
which  includes  a  zero  memory  nonlinearity.  We  analyze  the 
performance  of  this  nonlinear  receiver.  We  establish  the 
similarity  between  the  performance  of  the  nonlinear  receiver 
operating  in  combined  Cauchy,  Gaussian,  and  multiple  access  noise 
and  the  performance  of  a  linear  receiver  operating  in  combined 
Gaussian  and  multiple  access  noise. 


Another  contribution  of  the  thesis  is  to  extend  the  use  of 
the  characteristic  function  method  to  the  analysis  of  our  LF/CDMA 
signalling  scheme.  We  add  to  the  information  regarding  the 
shortcomings  of  the  Gaussian  approximation  for  CDMA  performance 
analysis . 

Another  contribution  of  the  thesis  is  the  application  of  a 
partial  sum  theorem  in  the  pursuit  of  a  minimax  sequence  design 
criterion.  Finally,  we  use  the  minimax  criterion  to  synthesize  2 
sequence  designs  for  our  LF/CDMA  scheme.  These  designs  provide  K 
sets  of  M  orthogonal  sequences.  Both  designs  have  "good" 
multiple  access,  skywave,  and  signal  acquisition  properties. 


CHAPTER  2 


SYSTEM  ANALYSIS  PART  I: 

RECEIVERS  WITH  NONLINEARITIES 

2.1  Introduction 

The  goal  of  Chapters  2  and  3  is  to  analyze  our  low  frequency 
CDMA  communication  system  and  present  equations  which  will 
describe  the  basic  system  tradeoffs.  To  this  end,  we  begin  this 
chapter  by  describing  our  LF/CDMA  communication  system  design 
model.  Next,  we  describe  LF  atmospheric  noise  and  describe  a 
noise  model,  which  is  appropriate  for  the  LF  channel.  We  then 
derive  a  locally  optimum  {small  signal)  Bayes  detector  for  M 
signals  and  approximate  the  performance  of  this  detector.  The 
most  interesting  feature  of  the  locally  optimum  Bayes  detector 
( LOBD )  is  a  zero  memory  nonlinearity  (ZNL)  and  the  performance 
analysis  shows  that  the  performance  of  the  overall  detector  is  a 
strong  function  of  Fisher's  information  number.  Consequently,  we 
conclude  Chapter  2  by  presenting  and  discussing  curves,  which 
describe  Fisher's  information  number.  In  Chapter  3,  we  continue 
our  system  analysis  by  considering  a  linear  analysis. 

2.2  Low  Frequency  CDMA  Commun icat ion  System  Design  Model 

In  this  section,  our  CDMA  design  model  is  presented  and 
discussed  in  detail.  As  shown  in  Figure  2.1,  every  user  has  a 
data  alphabet  of  size  M,  and  corresponding  to  each  element  of 
this  data  alphabet  is  a  length  N  signature  or  bandspreading 
sequence.  The  sequence  elements  are  rfc^ 


order  roots  of  unity; 


modulator 


2.1  LF/CDMA  Ccrrrnanication  System 


c'v 


1  N 


hence,  r  phases  are  used  to  modulate  the  carrier.  The  sequence 
set  or  code  for  the  user  is  denoted  by  ,  where 


C{ct)={ v(a)  v(a)  } 

’■vl  *  *  *  M  1 


(2.1) 


is  a  set  of  vectors. 


v(a)  ,] 
VA  ,N-1 J 


A=l,2, . . .M 


(2.2) 


The  sequence  corresponding  to  the  current  data  is  used  to  phase 
modulate  the  transmitted  carrier.  Therefore,  the  signal 
transmitted  by  the  transmitter  during  the  interval  t£[0,NT_) 


s  (v|a*  ,  t)  =Re{  I  r  (t-nT  )  expjui  t} 

A  n=0  A,n  c  c 


(2.3) 


r (t-nTc)  cos(wrt+9>  ) 
n=0  c  A'n 


where 


v  =expj6 
A ,n  A  ,n 


Ae{l,2...M} 
ne { 0 , 1 . . .N-l } 


(2.4) 


r(t)  is  the  chip  waveform  common  to  all  users.  It  is  of 
duration  Tc  and  is  the  envelope  of  the  modulated  carrier.  In  LF 
applications,  T(t)  need  not  be  a  constant  envelope  signal.  The 
energy  of  the  chip  waveform  is 


( t)  dt 


(2.5) 


•v.v.v 


(2.6) 


and,  consequently,  the  energy  of  the  signal  is 

NT 

ts=J  I  s  U)  (t)  i  2dt=Ne„  /2 
0 

The  carrier  frequency  shared  by  all  users  is  c oc  radians  per 
second.  Note  that  the  phase  of  the  carrier  during  the  nt'’1  chip 
interval  [nTc,  (n  +  1)  Tc)  is  controlled  by  the  n*-*1  sequence  element 
only.  With  LF  transmission,  no  additional  phase  offset  needs  to 
be  modelled.  In  other  words,  two  signals  will  be  in  phase  if 
their  relative  delay  is  zero  and  they  are  identically  modulated. 
The  bandwidth  expansion  factor  of  a  user's  sequence  set  is  the 
ratio  of  the  transmitted  bandwidth  to  the  information  rate 
(BW/R).  For  a  pulse  of  duration  Tc,  the  first  null  to  first  null 
transmitted  bandwidth  (BW)  is  approximately  2/Tc.  Therefore,  the 
bandwidth  expansion  factor  is  given  by 


BW/R=2N/log2M 


(2.7) 


Additionally,  the  transmitted  energy  per  information  bit  is 


£k=N€p/21og2M 


(2.8) 


In  the  channel,  each  signal  is  delayed  by  an  amount  t ^  and 
is  corrupted  by  atmospheric  noise  (nft(t)).  We  denote  the  vector 
of  signal  advances  by 


(2.9) 


In  this  thesis,  we  do  not  include  individual  attenuation  terms 
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for  each  signal,  but  these  could  easily  be  introduced.  Without 
loss  of  generality,  we  stud>  a  receiver  for  and  we  set  7(1) 

equal  to  zero,  since  we  assume  the  receiver  is  synchronized.  The 
received  signal  is  denoted  by 

r (t) =s (v^1 *  , t) +n (t)  te[0,NT  )  (2.10) 

where  the  data  word  sent  is  X  If  a  single  user  is 

signall ing , 

n(t)=nA(t)  (2.11) 

where  nA(t)  is  the  LF  noise  process  observed  at  the  input  of  the 

receiver.  If  more  than  one  user  is  transmitting, 

K 

n(t)=n  (t)  +  I  s(b(6),t+x(S))  t£  [  0  ,NT  )  (2.12) 

c 

S=2 

In  this  equation,  the  interfering  signals  are  shifted  in  time  by 
P  '  relative  to  the  reception  of  the  desired  signal. 
Consequently,  two  codewords  from  contribute  to  the 

modulation  of  the  signal  from  the  transmitter,  during  the 

interval  [  0 ,NTC) .  This  is  denoted  through  the  vector  b  )  (a  2N- 
tuple),  which  is  a  concatenation  of  the  first  and  second 
sequences  received  from  the  user  during  the  interval 

[0,NTc).  We  may  write 


This  notation  allows  us  to  write 


N-l 

s  (b(3)  ,t+T  (6)  )  =  I  Re{b]^  F  (t+T  (6)-nTc)expjwc(t+T  U)  )  }  (2.14) 

n=0 


where 


y t  <-(;)<(i+i)T 

c—  c 

Finally,  we  denote 


=expj  0 


(6) 

b,n 


(2.15) 


(2.16) 


In  the  next  section,  we  take  a  more  detailed  look  at  the  atmos¬ 
pheric  noise  term  nA(t). 


2.3  Atmospheric  Noise 

LF  atmospheric  noise  is  predominantly  generated  by  lightning 
strokes.  This  noise  has  much  greater  power  than  galactic  or 
receiver  noise  sources,  hence,  a  LF  noise  process  appears  as  a 
low  power  Gaussian  process  punctuated  by  high  amplitude  tran¬ 
sients.  Consequently,  the  atmospheric  noise  process  may  be  ex¬ 
pressed  as 


nA(t) =nG (t) +nc (t) 


(2.17) 


where  nG(t)  is  a  Gaussian  process,  and  nc(t)  is  an  impulsive 
noise  process.  The  Gaussian  component  is  due  to  the  combined 


contribution  of  noise  sources  which  are  either  weak  or  distant 


plus  thermal  noise  in  the  receiver.  In  our  receiver,  the  noise 
and  signals  pass  through  a  matched  filter  and  are  sampled.  The 
portion  of  this  sample  due  to  noise  is  denoted 


nA=nC+nG 


(2.18) 


where  V q  is  the  impulsive  component  and  where  ??q  has  a  Gaussian 
distribution.  The  "impulsiveness"  of  Vq  is  reflected  in  its 
voltage  deviation 

2  k 

vd=  (envelope (nA(t) )  ) }  /E{envelope (nA ( t) ) }  (2.19) 

The  voltage  deviation  of  band  pass  atmospheric  noise  is  plotted 
in  Figure  2.2  as  a  function  of  center  frequency  and  bandwidth. 
For  wide  bandwidths,  the  deviation  is  large,  because  a  small 
number  of  atmospheric  noise  events  will  dominate  the  output  of 
the  filter.  On  the  other  hand,  for  narrow  bandwidths,  the  devia¬ 
tion  decreases  and  approaches  the  value  for  a  Gaussian  process. 
This  is  because  the  filter's  time  constant  is  large  compared  to 
the  interarrival  time  of  the  impulsive  events  and  the 
contribution  of  any  one  impulse  cannot  dominate.  In  Figure  2.2, 
the  deviation  curves  for  two  sets  of  bandwidth  are  shown.  The 
upper  curves  correspond  to  bandwidths  that  would  be  typical  of 
our  receiver's  front  end  and  the  lower  curves  correspond  to  the 
post  processing  bandwidth  when  sequence  lengths  near  1000  are 
used.  A  more  thorough  investigation  of  the  noise  statistics  is 
required  to  be  able  to  predict  signal  detector  performance. 


VOLTAGE  DEVIATION 


Atmospheric  noise  has  been  characterized  through  its  first 
order  probability  density  function  (pdf)  in  a  number  of  ways  in 
the  recent  literature.  The  noise  pdf  has  been  treated  as  a 
generalized  Cauchy  distribution  by  Hall  (1966)  and  Feldman 
(1972).  Also,  a  log  normal  pdf  has  been  used  by  Omura  and  Shaft 
(1971).  However,  we  will  review  the  results  of  Middleton  (1974, 
1976,  1977,  1979a,  1979b)  which  are  attractive,  because  the 
mathematical  form  of  his  results  does  not  change  with  changing 
physical  conditions.  However,  the  parameters  of  the  model  are 
explicit  functions  of  the  underlying  physical  mechanisms 
(physical  source  distributions,  noise  source  waveforms, 
propagation  types,  etc.).  This  means  that  the  important 
parameters  of  the  model  may  be  readily  measured.  It  also  means 
that  the  connection  between  analytic  simplifications  and  the  real 
world  may  be  established.  Middleton's  noise  model  is  limited  to 
the  description  of  noise  after  it  has  been  processed  by  a 
narrowband  filter.  (In  this  case,  this  means  that  the  receiver 
front  end  bandwidth  must  be  less  than  one  fifth  of  the  center 
frequency.)  The  other  model  limitation  is  that  the  noise  must 
only  be  statistically  related  to  the  receiver.  (Deterministic  or 
completely  known  interference  must  be  treated  in  another 
fashion . ) 

Three  classes  of  noise  are  defined  by  Middleton.  Class  A 
noise  is  where  the  unfiltered  noise  has  a  power  spectral  density 
(psd)  roughly  equal  to  or  less  than  the  receiver  front  end  band¬ 
width.  Class  B  noise  is  where  the  unfiltered  noise  process  has  a 
psd  substantially  wider  than  the  receiver  front  end,  and  class  C 


noise  is  a  sum  of  a  class  A  and  a  class  B  noise  process.  The 
atmospheric  noise,  which  interests  us,  is  clearly  a  class  B  type 
noise  process.  For  all  three  noise  types,  Middleton  has  devel¬ 
oped  a  characteristic  function,  envelope  exceedance  probability, 
envelope  pdf  and  envelope  moment  expression.  All  three  of  these 
expressions  come  in  pairs  for  class  B  noise,  because  one  expres¬ 
sion  is  required  to  describe  the  small  envelope  noise  behavior 
and  another  is  required  for  the  large  envelope  behavior.  For 
Middleton's  full  class  B  model,  six  parameters  are  involved; 
however  Spaulding  (1982)  has  shown  that  for  signal  detection  work 
in  atmospheric  noise,  a  simplified  version  of  the  model  is  quite 
satisfactory.  In  other  words,  we  will  only  consider  the  small 
envelope  forms  of  the  class  B  model.  The  small  envelope  pdf  for 
class  B  noise  is: 

exp(-z2/k1) 

p  (2)  = -  .  (2.20) 

'  C 

/  i  >  m 
00  (-1) 

I  -  A™  T  ( (ms+l)/2)  ,F.  (— ms/2,?5,z2/k1  ) 

m=0  ml  s  1  *  1  1 

where  is  the  Gamma  function  and  ^F^  is  a  confluent 

hypergeometric  function  described  by  Middleton  (1976).  The  pdf 
has  three  parameters:  s,  As  and  k  ,  where  s  and  As  are  intimately 
linked  to  the  physics  of  the  noise  process,  and  k-^  is  a 
normalization  parameter  such  that  the  variance  of  the  process 
envelope  is  unity.  The  variable  s  is  called  the  spatial  density 
propagation  parameter  and  describes  the  spatial  distribution  of 
the  noise  sources  as  well  as  the  noise  source  to  receiver 


propagation  law.  This  parameter  is  limited  to  the  values  0<s<2 
in  Middleton's  analysis,  where  he  states  that  these  values  cover 
most  of  the  practical  cases.  As  is  the  effective  impulsive  index 
of  the  noise  and  is  a  function  of  the  impulsive  index  of  the 
noise  (AB)  as  well  as  s.  The  impulsive  index  (AB)  is  equal  to 
the  mean  number  of  interference  events  per  second  times  the  mean 
duration  of  an  event,  where  the  duration  is  measured  after  the 
front  end  filter.  As  such,  AB  is  a  measure  of  the  temporal 
overlap  or  density  of  the  noise.  As  the  mean  number  of 
interfering  events  per  second  increases,  Afi  and  As  increase  and 
the  noise  pdf  resembles  a  Gaussian  pdf  as  might  be  expected, 
because  of  the  central  limit  theorem. 

For  s  -  1,  Middleton's  class  B  pdf  may  be  approximately 
simplified  (Middleton  (1976))  to  a  Cauchy  pdf 


2a, 

Pn  2  2 

nc  ir(4z  +k1A£) 


(2.21) 


with  effective  impulsive  index  Aj_,  and  normalization  k -j_ .  A 
revealing  expression  for  A-^  is 


Al=4v«;/va 


(2.22) 


where  v^  is  the  voltage  deviation  given  by  Equation  (2.19)  and  AB 
is  the  impulsive  index  described  above. 


In  conclusion,  the  pdf  of  Equation  (2.21)  is  not  general 
enough  to  describe  all  varieties  of  atmospheric  noise,  but  it  is 
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a  conservative  choice  for  the  results  we  pursue  here  and 
consequently  we  will  use  it  in  the  remainder  of  this  chapter.  We 
now  turn  our  attention  to  the  design  and  analysis  of  locally 
optimum  signal  detectors. 


2.4  A  Locally  Optimum  Bayes  1  Detector  of  M  Signals 

As  mentioned  in  Section  2.3,  low  frequency  atmospheric  noise 
has  a  distinctly  non  Gaussian  nature.  Specifically,  samples  of 
atmospheric  noise  after  front  end  filtering  may  be  denoted 


nA=nC+riG 


(2.23) 


where  V c  is  the  predominant  component  and  has  a  Cauchy  first 
order  pdf.  The  performance  of  a  signal  detector  can  be  a  strong 
function  of  the  "match"  between  the  assumed  noise  model  and  the 
actual  noise,  consequently,  we  wish  to  investigate  which  signal 
detection  scheme  is  optimum  for  our  noise  model.  Our  optimized 
detector  should  give  better  performance  than  the  optimum  Gaussian 
noise  detector  operating  in  Gaussian  noise,  because  Gaussian 
noise  impairs  detector  performance  more  than  any  other  noise  type 
(for  a  given  noise  variance).  The  analysis  of  signal  detection 
in  non-Gauss i an  noise  is  a  widely  studied  subject  and  some  well 
known  results  are  published  in  Miller  and  Thomas  (1972), 
Spaulding  and  Middleton  (1977),  Spaulding  (1982),  Omura  and  Shaft 
(1971),  and  Lu  and  Eisenstein  (1981). 

In  Spaulding  and  Middleton  (1977),  the  optimum  coherent 


detection  of  binary  sianals  in  class  A  impulsive  noise  is 


discussed  in  three  main  parts.  First,  the  optimum  Bayes' 
decision  strategy  is  derived  for  class  A  noise  and  an  average 
probability  of  error  bound  derived.  Second,  the  performance  of 
sub-optimum  matched  filter  detectors  in  class  A  noise  is 
calculated.  In  the  third  section,  the  optimum  strategy  derived 
in  the  first  part  is  simplified  through  the  use  of  a  small  signal 
assumption.  This  simplification  is  important,  because  the  rather 
complex  globally  optimum  test  is  reduced  to  a  matched  filter 
preceded  by  a  zero  memory  nonl inear i ty  as  derived  earlier  by 
Miller  and  Thomas  (1972).  This  theoretically  derived  ZNL/matched 
filter  structure  is  interesting,  because  it  resembles  earlier 
intuitively  derived  receiver  designs.  The  intuitive  receiver 
designs  for  impulsive  noise  were  based  on  the  observation  that 
limiting  the  amplitude  of  the  received  signal  will  improve  signal 
detectability,  because  the  majority  of  the  noise  power  is  con¬ 
tained  in  high  amplitude,  short  duration  transients. 

The  small  signal  assumption  can  be  used  to  derive  locally 
optimum  detectors  for  signals  imbedded  in  any  variety  of  noise 
process  and  it  can  be  used  to  approximately  analyze  the 
performance  of  these  detectors.  Indeed,  this  assumption  has  been 
used  to  derive  and  analyze  detectors  for  signals  corrupted  by 
class  B  noise  processes  (Spaulding,  1982).  In  this  section,  we 
will  use  it  for  our  model  where  the  signal  is  degraded  by 
additive  impulsive  noise,  additive  Gaussian  noise  and  multiple 
access  interference.  As  such,  the  small  signal  assumption  is 
well  suited  for  use  with  our  model  where  the  desired  signal  is 
generally  small  compared  to  the  sum  of  the  degrading  signals. 


The  locally  optimum  (small  signal)  Bayes  detector  for  our 
model  is  now  derived  and  its  performance  approx imated.  A  sampled 
data  receiver  is  assumed  where  the  sampler  is  preceded  by  a 
correlation  receiver  which  also  serves  to  reduce  out  of  band 
interference  and  reduce  dynamic  range  problems  in  analog  to 
digital  conversion.  The  local  signal  in  the  receiver  is  the 
carrier  (with  zero  phase  shift)  modulated  by  the  chip  waveform 
(i.e.,  r(t)coswct).  This  situation  is  shown  in  Figure  2.3, 
where  the  received  signal  is 

r(t)=s^  (t)+n(t)  (2.24) 


The  sampled  output  of  the  receiver  at  time  t  =  nTc  is 

nTc  nTc 

R  =  f  n(t)  nt)cosw  t  dt  +  f  s  (1)  (t)  T  (t)  cosu>  t  dt  (2.25) 

n  j  c  *  c 

(n-l)Tc  (n-l)Tc 

=n+  (<or/2)  cos6,(1^ 
n  i  A ,  n 

where  V n  is  the  result  of  the  atmospheric  noise  plus  multiple 
access  interference.  The  set  j(Rn)|  are  the  observations  from 

which  our  detector  must  make  an  optimal  decision. 

The  Bayes  decision  criterion  leads  to  a  likelihood  ratio 
test  (LRT)  ,  and  when  the  signals  are  equally  likely  and  symmetric 
cost  assignments  are  made,  the  LRT  is 


choose  A 


P 


,  (x  !  A  )  >  p_  , (x  !  u) 


(2.26) 


We  denote 


R= ( Rq  ...  Hn_1>  (2.27) 

n=(nQ  ...  (2.28) 

S  (1)  =  (J5£r)  (cosB^q - cosB ^ )  (2.29) 

as  the  observation,  noise  and  signal  vectors,  respectively.  We 

now  observe  that 

pR|X(x|X)=pn(x-sj[1))  (2.30) 


which  we  approximate  near  .  5ep  =  0 ,  for  threshold  signals  as 
follows 


p  (x) 


(2.31) 


N-l 

-  I 

n=0 


9x 


(JjEp)  COS0 


(1) 
A ,  n 


Assuming  independent  noise  samples,  along  with  manipulating 
Equations  (2.26)  and  (2.31)  yields  an  equivalent  LRT 

(2.32) 

V  U 

where  a  zero 
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Z 


(1) 


n=0  dR 


N-l 
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n 

d 


log  p  (R  )  c osB'*'  < 
hn  n  A  ,n 


(1) 


n=0  dR_ 


n 


log  p  (R  )  cos9 

nn  n  h  ,  n 


A  receiver  with  this  LOBD  is  shown  in  Figure  2.3, 


Figure  2.3  A  Locally  Optimum  Detector  of  M  Signals 
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memory  nonlinearity  maps  the  samples  Rn  into  Rn.  The  equation  for 
the  zero  memory  nonlinearity  (ZNL)  is 


/\ 


d 

dR 


n 


log  pv(R  ) 

In  “ 


(2.33) 


The  Rn  are  then  processed  by  M  discrete  time  correlators  each 
matched  to  one  of  the  M  signals.  The  receiver  decides  which  data 
word  was  sent  by  choosing  the  data  word  that  corresponds  to  the 
correlator  with  the  smallest  output. 

We  now  approximately  analyze  the  performance  of  the  LOBD,  by 
defining  the  signal  to  noise  ratio  as  follows. 


SNR ( y ,  A )  = 


E  {*1  |  X  }  -E  K  J  A  } 
U  A 


(2.34) 


(Var  (/Cj  A)+Var(/tx  |  A  )p 
If  the  J7n  are  Gaussian  (or  the  and  are  uncorrelated  and 
Gaussian),  then  this  definition  is  especially  appealing.  For  the 
Gaussian  case,  the  probability  of  codeword  error  can  be  union 
bounded  as  follows 


M 

Pr(e)  <  (1/M)  I  I  Q (SNR)  (2.35) 

A=1  u/A 


where 

1  00 

Q(z)  =  ~ZI-  /  exp(-y2/2)  dy 


(2.36) 
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The  moments  of  can  be  found  as  functions  of  the  moments 


of  &ncos#i^ 

n  A,n 


E { R  cos0 ,  | X  sent }  = 

n  A ,  n 


(2.37) 


E(1){  Cr  cos 0,(1I  pn  (R  -S.  n)  dR  . 
c  J  n  A,n  n  A,n  n 


In  Equation  (2.37),  Ec ^  denotes  expectation  with  respect  to  the 
phase  of  the  nth  chip  (0^).  In  what  follows,  this  phase  is 
approximated  as  a  uniform  random  variable  on  [0,27T  ),  because  the 
number  of  phases  (r)  in  our  LF/CDMA  system  is  typically  quite 
large.  This  approximation  is  further  discussed  in  Section  3.3. 
Employing  the  small  signal  approximation  for  P77n(Rn)  yields 


E{R  cosO,  | A  sent}  = 
n  A  ,n 


-E  ( 1 }  {  ( cos0  ! )  2  L  Jser} 
c  A ,  n  i 


(2.38) 


where 


■dly 1x1)2 

1  lx) 

n 


(2.39) 


The  quantity  L  is  known  as  Fisher's  information  number 
Add i t i ona 1  ly 


E { R  COS0,  | A  not  sent}  = 
n  A  ,  n 


0 


(2.40) 


To  find  the  variance,  the  following  moments  are  considered 


~  111  1 

E { ( R  cos9,  ')  X  sent }=  (2.41) 

n  a  f  n 

E(1)  f  (R  )  2  (COS0,(1)  )  2  pn  (R  -S,  )  dR  =E{1)  (coseJ1^)2L  = 

c  J  n  X,n  nA,n  nc  A,n 

-oo  n 

E{ (R  cos0 | ^ ) 2 | X  not  sent}  =  hL  (2.42) 

n  a  f  n 

Hence , 

Var  (R^osQ^1^  |  X  sent)  =  (%L)  -  (hLe  p)  2  (2.43) 

~  JsL 

because  of  our  small  signal  approximation. 

Var(Rncos0^  ^|X  not  sent)  =  hL  (2.44) 

The  moments  of  and  are  now  given  by 

E{/^|A  sent}  =  -!sNLer  45^ 

E{<’u|X  sent}  =  0  (2.46) 

VarK\J\  sent }  -  Varf.f^  j  X  sent}  =  NL/2  (2.47) 

Consequently,  the  signal  to  noise  ratio  is 


If  the  interfering  noise  process  n(t)  was  a  Gaussian  process, 


then  the  samples  would  have  Gaussian  distribution  and  the 
solution  of  Equation  (2.36)  would  yield 


L  =  (Var (n  ) )  1  =  a~2 
n  G 


(2.49) 


Additionally,  the  receiver  would  be  completely  linear,  because 


R  =  -R 
n  n 


(2.50) 


Finally,  the  signal  to  noise  ratio  would  become 


SNR  =  (Np  !ser/oG 


(2.51) 


Equation  (2.51)  causes  us  to  define  the  processing  gain  of  the 
zero  memory  nonlinearity  as  follows 

ZNLPG  =  L  var (n  )  (2. 52) 

n 


Since  our  input  noise  has  unbounded  variance  the  ZNL  "processing 
gain"  is  also  unbounded.  Clearly  a  ZNL  of  some  sort  will  always 
be  required  in  CDMA  systems  suffering  from  impulsive  noise.  in 
the  next  section,  we  will  calculate  L  for  our  noise  model. 


2.5  Performance  in  Combined  Atmospher ic  No i se  and  Multiple 


Access  Interference 

In  this  section,  we  will  calculate  L  as  a  function  of  the 
parameters  in  our  noise  model.  From  Equation  (2.39),  L  is 


,  d  ,  .  .  2 
( j —  p  (x)  ) 

”  dx  *nn 

L  =  J  -  dx  (2.53) 

P  (x) 
nn 

where  ?7n  is  the  n*-1"1  noise  sample.  This  expression  is  difficult 
to  evaluate  directly,  so  we  integrate  by  parts  and  take  limits  to 
obta  i  n 


L 


log  p  (x)  r  p  (x)  dx 
nn  dx^  nn 


(2.54) 


Symmetry  considerations  provide 


L 


p  (x) dx 
nn 


(2.55) 


The  characteristic  function  of  the  random  variable  77n  may  be 
defined  as  the  Fourier  transform  of  the  first  order  probability 
distribution  function  and  is  denoted  <£>^(u>).  Characteristic 
functions  are  one  way  to  obtain  an  expression  for  L  that  can  be 
readily  evaluated  via  computer.  Consequently,  we  observe  that 


(x)  =  (1/tt)  f  p  (u>)  cosfjx  do 


n 


(2.56) 


.  2  pn 
dx  n 


=  (-1/7T)  /b  (o)o' 


co  sox  d'. 


(x) 


(2.57) 
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Each  noise  sample  is  a  sum  as  follows 


■  "c  *  Bh  n 


(2.58) 


where  77^^  is  the  portion  of  the  noise  due  to  the 
competing  user.  Since  each  term  on  the  right  hand  side  of 
Equation  (2.58)  is  independent,  the  character istic  function 
0^(o>)  is  a  product 


/ o \  k— 1 

<p  (oj )  =<p„  (w)  <pr  (d)  (4>  ;  (w) ) 

n  C. 


(2.59) 


If  the  individual  characteristic  functions  are  known,  the  ZNL 
processing  gain  may  be  computed. 

As  discussed  in  Section  2.3,  the  pdf  of  Vq  is  being  modelled 
as  a  Cauchy  pdf.  Consequently, 


4>c((jj)  =  exp  (- lwA^/2  )  ) 


(2.60) 


The  Gaussian  portion  of  the  noise  has  pdf 


p  (z)  =  — — 
nG  «^2tt  o. 


exp (~z2/2°^) 


(2.61) 


and  characteristic  function 


(t>G('jj)  =  exp(- (^aQ)  V2) 


(2.62) 
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Finally,  the  multiple  access  interference  introduced  by  each 
competing  user  is 


n(8)  -  cose*8' 


(2.63) 


where  O'  may  be  modelled  as  a  uniform  random  variable  on 
( 0 , 2  7T )  as  discussed  in  Section  3.3.  Eachf/^)  has  the  following 


distribution 


pn<6) 


(z)  = 


—  <1- <  2z/er )  2)  _>s 


-HZT<Z<HCy 


otherwise 


(2.64) 


The  variance  of  this  noise 


and  characteristic  function  is 


<t> (w)  =  JQ(w/2a^;  ) 


(2.65) 


(2.66) 


where  J0  is  the  Bessel  function  of  order  zero. 

The  above  equations  have  been  used  to  plot  a  normalized 
Fisher's  information  number  defined  as  follows 


L  =  (var(nG)  +  (K-l)var(n  ))  L 


(2 . 6‘ 
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In  Figures  2.4  through  2.8,  L  is  plotted  versus  the  variance  of 
the  Gaussian  noise  component  and  with  the  number  of  users  (K)  as 
a  parameter.  Figure  2.4  is  for  A  j.  =  .  0  1  and  v  a  r  ( i)  O  )  )  =  #  g  g  i  < 
Figures  2.5,  2.6,  2.7,  and  2.8  are  for  A^  =  .l  and  var  (7^0)  )  =0.001 , 
0.01,  0.1,  and  1.0  respectively. 

The  above  equations  have  also  been  used  to  plot  the  form  of 
the  ZNL  given  by 

R  =  -  -  log  p  (R  )  (2.68) 

n  dR  nn  n 

n 

In  Figures  2.9  and  2.10,  the  ZNL  is  plotted  for  K  =  ll, 
var  iance  (  =.l,  and  A^  =  .01  and  A^  =  .l  respectively. 

We  now  discuss  the  figures.  The  most  important  feature  of 
Figures  2.4  through  2.8  is  that  L  approaches  1  with  increasing  K 
or  Gaussian  noise  variance.  In  other  words,  as  multiple  access 
noise  power  ££  Gaussian  noise  power  increases,  Fisher's 
information  number  approaches  their  combined  variance.  This 
means  that  the  LOBD  performance  is  close  to  the  performance  of  a 
linear  receiver  without  the  impulsive  noise.  To  see  why, 
consider  the  ZNLs  plotted  in  Figures  2.9  and  2.10.  All  these 
optimal  ZNLs  are  linear  for  small  signals  and  strongly  suppress 
signals  above  a  certain  threshold.  The  threshold  value  is 
approximately  3.5  times  the  standard  deviation  of  the  combined 
multiple  access  and  Gaussian  noise.  As  the  multiple  access  or 
Gaussian  noise  power  increases,  the  receiver  spends  more  of  its 
time  operating  in  the  linear  region.  At  the  same  time,  the 
impulsive  events  that  do  occur  are  suopressed  by  the 


Figure  2.6  Normalized  Fisher's  Information  Number  versus 
Variance  of  Gaussian  Noise  Component 


Normalized  Fisher's  Information  Number  versus 
Variance  of  Gaussian  Noise  Component 
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Figure  2.9  Optimal  ZNL  for  Combined  Cauchy,  Gaussian 
and  Multiple  Access  Noise 


nonlinearity.  From  the  above  observations,  we  draw  two  major 
conclusions . 

First  of  all,  every  LF/CDMA  receiver  should  employ  some  form 
of  nonlinearity  to  suppress  the  effect  of  the  impulsive  noise. 
This  nonlinearity  need  not  be  as  complicated  as  the  ones  given  in 
Figures  2.9  or  2.10.  Certainly,  a  simple  clipper  or  hole  puncher 
will  work  nearly  as  well.  These  simpler  ZNLs  could  be  designed 
to  suppress  inputs  greater  than  3.5  times  the  standard  deviation 
of  the  combined  multiple  access  and  Gaussian  noise. 
Alternatively,  the  threshold  could  be  designed  to  adapt  such  that 
a  certain  percentage  of  the  input  samples  where  suppressed 
(Feldman  (1972)). 

Our  second  major  conclusion  is  that  the  LOBD  operating  in 
combined  multiple  access,  Gaussian  and  impulsive  noise  has 
performance  close  to  a  linear  receiver  operating  in  the  same 
environment  without  the  impulsive  noise.  The  Gaussian  noise 
power  in  the  linear  analysis  could  be  adjusted  to  include  the 
integrated  effect  of  a  large  number  of  clipped  impulsive  events. 


CHAPTER  3 


SYSTEM  ANALYSIS  PART  II: 
LINEAR  RECEIVERS 


3.1  Introduction 

In  this  chapter,  we  consider  linear  receivers  operating  in 
the  presence  of  additive  white  Gaussian  noise  (AWGN) .  In  the 
first  section,  we  will  only  consider  a  single  user  system,  but  in 
Section  3.4,  we  will  consider  multiple  users.  Between  Sections 

3.2  and  3.4,  we  will  discuss  the  random  sequence  analysis  ideas, 
which  are  important  to  our  multiuser  analysis. 

3.2  Single  User  Systems 

Our  linear  receiver  is  shown  in  Figure  3.1  and  it  simply 
consists  of  M  correlators  and  a  "largest  of"  decision  device. 
Each  correlator  is  matched  to  one  of  the  signals  from  the  code 
and  if  K  equals  one 

INe  r 

— - —  +  I...  if  A  sent  (3.1) 

2  A*  A 

t  ,  if  A  not  sent 

A,  A 

IA^  is  the  inner  product  of  s(v^^,t)  and  the  AWGN  process, 
which  has  power  spectral  density  equal  to  NQ/2. 

The  performance  of  our  linear  receiver  (in  Figure  3.1)  for 
single  users  and  AWGN  is  analyzed  in  many  communication  texts 
such  as  VanTrees  (1968)  and  Wozencraft  (196  5).  The  union  bound 
can  be  used  to  provide  the  following  well-known  probability  of 


codeword  error  bounds 


Q(/es/Nq)  <Pr  (c)  <  (M-l)Q(/es/No) 


(3.2) 


where  the  codewords  are  orthogonal  and  Q  is  given  in  Equation 
(2.41).  For  comparing  sequence  sets  of  different  sizes,  a  bound 
for  Pr (e)  as  a  function  of  the  energy  per  information  bit  eb  is 
more  useful.  Such  an  expression  is  made  possible  by  using 

£b=es/l0g2M  (3.3) 

Now, 

Q(/eblog2M/NQ)  £Pr  (e)<_  (M-l)  Q  (/eblog2M/NQ)  (3.4) 


Throughout  this  thesis,  we  will  assume  each  user's  code  is 
orthogonal,  and  we  will  now  briefly  describe  a  construction  for 
these  codes,  which  is  suitable  for  a  single  user  system.  An 
orthogonal  code  may  be  constructed  as  follows 


C 


(0) 


M 


(v, 

A 


*=1 


(3.5) 


V  i  ~  (  V  p  ...  V. 

A  A  t  U  /  #  1 

v,  -exp ( j  2 " * n/N ) 

a  ,  n 


(3.6) 


(3.7) 


This  code  may  be  represented  by  a  matrix  D,  whose  rows  are  the 
code  words 


D= 


V0,0  V0,l  *  *  V0,N-1 


l_VN-l,0  -  VN-1,N-1J 


(3.8) 


The  matrix  D  is  the  discrete  Fourier  transform  matrix  and  the 
corresponding  orthogonal  codes  have  the  following  properties 

r=N  (3.9) 

M=N 


Orthogonal  codes  (constructed  in  any  way)  are  characterized  by 
their  probability  of  error  bound  versus  £b/N0  ratio  in  Figure 
3.2,  where  sequence  length  (or  equivalently,  code  size)  is  a 
parameter . 

3.3  Random  Sequence  Analysis 

In  the  last  section,  we  analyzed  the  performance  of  our 
linear  receiver  for  the  case  of  one  user,  and  in  the  next 
section,  we  will  expand  our  analysis  to  the  multiuser  case.  For 
part  of  our  multiuser  analysis,  we  will  use  random  sequence 
analysis.  Hence,  we  use  this  section  to  continue  the  discussion 
of  random  sequence  analysis,  which  we  began  in  Section  1.2.  As 
mentioned,  taking  the  expected  value  of  system  performance 
parameters  with  respect  to  the  random  sequences  makes  it  possible 
to  get  useful  approximate  results  without  designing  any  actual 


codes.  The  role  of  random  analysis  in  the  analysis  of  CDMA 
communication  systems  is  explained  in  Roefs  and  Pursley  (1976) 
and  Pursley  (1977). 

For  some  of  our  multiuser  analysis,  each  user's  code  will  be 
treated  as  a  random  coset  of  an  orthogonal  code.  In  other  words, 
each  codeword  in  consists  of  the  element  by  element  product 
of  a  codeword  from  an  orthogonal  code  and  a  random  vector  (or 
coset  leader).  Mathematically,  if 

C(0,-{x<0))M  <3.10) 
A  A=1 


is  an  orthogonal  code,  then  every  codeword  in  may  be 

expressed  as 
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(a) 

X 


,x 


(O) 

X  ,N-1  * 


c 


(a)  , 

N  — 1  ’ 


(3.11) 


where  the  are  independent  identically  distributed  random 


variables.  Additionally, 


n 


(3.12) 


(ai 

For  our  work,  the  9C/n  will  be  uni  f  ormly  d  i  str  ibuted  from  0  to 
27 T,  because  r  is  so  large  that  a  uniform  distribution  is  a  good 
approximation.  Expectation  with  respect  to  the  random  vector 
c^  will  be  denoted  eJ°^  . 

Random  sequence  analysis  will  be  used  in  the  next  section  to 
approximate  the  probability  of  error  for  our  LF/CDMA  model  as  a 


function  of  the  number  of  users. 


m 


3.4  Multiuser  Analysis 


3.4.1  Introduction 


In  this  subsection,  we  develop  a  probability  of  error 
estimate  and  bound  for  the  multiuser  system  model.  The  estimate 
is  based  on  the  assumption  that  the  interference  introduced  by 
the  multiple  users  has  a  Gaussian  distribution.  Consequently,  it 
only  requires  the  calculation  of  the  variance  of  the  multiple 
access  interference.  Second,  we  bound  the  system  probability  of 
error  via  the  characteristic  function  method.  This  technique  is 
pioneered  for  CDMA  systems  in  Geraniotis  and  Pursley  (1982)  and 
Geraniotis  (1983)  and  we  obtain  some  interesting  extensions. 

For  our  multiuser  analysis,  we  will  continue  to  use  the 
receiver  of  Figure  3.1,  but  now  the  output  of  each  correlator  is 


given  by 


r  c 

— 2~  *  1A,\*  '\ 


X  sent 


X  not  sent 


(3.13) 


where  is  the  multiple  access  interference.  The  multiple 
access  interference  is  the  sum  of  interference  introduced  by 


each  of  the  competing  users: 


where 


1\(i)=f0  C  sCv^15  ,t)  s  (b(B)  ,t+x  (B)  )dt 


(3.15) 


To  rewrite  Equation  3.15  and  to  make  analysis  simpler  we 
introduce  some  new  functions.  First  of  all,  the  correlation 
between  the  signature  sequence  x  and  the  concatenation  of  the 
sequences  u  and  v  is 


H 


x,u,v 


.  i,-l 

U>«  I 

n~  X 


X  ,  u  + 
n-1  n 


1-1 

I 

n=0 


x  i  |  v 
n+N-J^  n 


(3.16) 


This  function  may  be  related  to  the  aperiodic  cross-correlation 
functions  described  in  Sarwate  and  Pursley  (1980): 


Hx,u,v<*'-Cx,u(1>+Cv,x(N-*> 


(3.17) 


Also  note  that 


Hx,u,u'S»=9x,u«» 


(3.18) 


where  0XfU(O  periodic  cross-correlation  function  which  is 

also  extensively  discussed  in  Sarwate  and  Pursley  (1980).  In 
addition,  the  partial  autocorrelation  functions  of  the  chip 
waveform  (T(t))  are  given  by  (Pursley,  1982) 

T 

ftr(s)=jC  r  (t)  r  (t-s)dt  0<s<Tc  (3.19) 

s 


and 

s 

Rr(s)  =  f  r(t)T(t+T  -s)  dt  0<s<T 
I  c  —  —  c 


(3.20) 


For  our  preferred  example  chip  waveform  (the  sine  pulse),  these 
functions  may  be  developed  to  yield  (Pursley,  1982) 


1L  (s)  =_£_  (  (T  -s)  cos  (tts/T  )  +  (T  /tt )  sin  (tts/T  )) 
1  ip  c  c  c  c 

c 


(3.21) 


Rp(s)= — (-scos(irs/T  )  +  (T  /tt)  sin  (tts/T  )) 


(3.22) 


Note  that 


ftr(0)-Rr(Tc)»er 


(3.23) 


Rr (Tc)=Rp(0)=0 


(3.24) 


With  the  correlation  functions  defined  above  and  equation  2.15, 
we  may  develop  Equation  3.15  as  follows 


I  (8)=JsR(t  (B)  )  (cosio  x  m  Re{H  ({)  } 

A  C  X  t  U  /  V 


(3) 


(3.25) 


-sinw  t (6)  Im{H  „  (5L)  )) 

C  X  ,  U  ,  V 


+  hR(T  ^  )  (COSO)  T  ^  Re{H  11  v(X+1)} 

C  A  f  u  f  V 


(S) 


-sinw  x  (6)  Im{H  V(&+1)  >) 

C  X  f  u  /  v 
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or 


I^6)=i5ft(T  (6)  )  Reiexpju^x  (6)  (C*  ,  u  ($  )  +Cv> x  (N-\)  )  } 


(3.26) 


+j5R(t  (e)  )  Re(expju)  x  l£5^  (c"  (l+l)+C  (N-|-l))} 

C  A  ^  U  V  /  A 


(B) 


where 


b(8)=(u(6)|v(6)) 


XT  ,<T  (e)  <(S+1)T 


(3.27) 


(3.28) 


As  shown  above,  the  interference  introduced  by  each  of  the 
competing  users  is  a  function  of  the  sequences  sent  by  that  user 
during  the  integration  time  [0,NTc)  as  well  as  the  relative 
delay.  In  this  thesis,  both  the  choice  of  codewords  and  the 
relative  delay  are  modelled  as  random.  Specifically,  the  u^ 
and  v  are  independent  and  equally  likely  to  be  any  sequences 
in  C<0>.  Additionally,  u^)  and  v^)  are  independent  of  u^^ 
and  vM  for  all  /3  not  equal  to  7.  Expectation  with  respect  to 
the  /3th  transmitter's  random  choice  of  data  is  denoted 


E*6)  (G(b(B)  )  }=(1/M2)  I 


utC(S>  vsC(6>  G<b 


(B) 


(3.29) 


Expectation  with  respect  to  the  choice  of  data  by  all  the 
competing  users  is  denoted  and  the  corresponding  variance  is 
denoted  var^. 
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The  relative  signal  delay  is  treated  as  a  random  variable, 
which  is  uniformly  distributed  over  the  interval  [0,NTc).  We 
denote  expectation  over  this  random  variable  as 


NT 


{G(T (B) ) }-(l/NTc) J  G(T(6))dT 


(8) 


(3.30) 


0 


Expectation,  with  respect  to  all  the  competing  user  signal 
delays,  is  denoted  E  j  and  the  corresponding  variance  is  varT  . 
The  average  probability  of  error  for  our  receiver  is 

M 


Pr(e)  =  (l/M)  X  Pr  (error  I  sent) 

X=1  A 


(3.31) 


Employing  the  union  bound  yields 


M 


Pr(e)<(l/M)  I  Z  Pr  (/?>(fj  vx 
X=1  y^A 


v.(1)  sent) 


(3.32) 


The  union  bound  upper  bounds  the  probability  of  the  union  of 
some  number  of  events  by  ignoring  the  probability  of  the  various 
joint  events.  Consequently,  it  is  a  tight  bound  when  the  sum  of 
probabilities  is  small  (less  than  10“^)  and  it  is  not  so  tight 
when  the  probability  of  error  is  large.  Substituting  expressions 


for  and  yields 


M 


?r(e)<(l/M)  I 


X  =  1  U^X 


Pr ( I  , -I.  . 
A, y  A, A 


Ner 

+  'cl  x>  I  A  sent) 


(3.33) 


M  2 
= (1/M)  I  I  Pr ( — 

A=1  y^A  Ne, 


(I.  -I.  ,  + )  >l|  X  sent)  (3.34) 

A ,  u  A ,  A  u  a 


Lam 
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The  right  hand  side  of  this  inequality  is  developed  by  two 
different  techniques  in  the  next  two  subsections. 

3.4.2  Gaussian  approximation 

The  inequality  of  3.33  becomes  an  approximation  when  the 
following  observations  and  assumptions  are  applied.  First, 
observe  that  *s  a  Gaussian  random  variable  with  zero 

mean  and  variance  N£pNQ/2  and  it  is  independent  of  the  multiple 
access  interference.  Second,  the  random  variable  (tij ^as 
zero  mean  and  its  distribution  is  assumed  to  be  Gaussian,  which 
allows  us  to  write 


M 

• 

"2No  4 

Pr(e)^(l/M)  I 

I  Q 

+  2  2  varb,T^'i^u 

A=1 

Ner  N  Ej/  D,T  U 

- 

As  shown  in  Equation  (3.22),  the  multiple  access  interference 
contributes  to  the  system  performance  through  its  own  second 
moment.  Due  to  the  independence  of  the  interference  introduced 
by  the  different  users,  we  may  write 


v*rb  T  (n£  -^A)=  I  E}!6)E^B)  (I  -I  ) 
b,T  p  A  g=2  b  x  v  X 


(3.36) 


/ 


If  the  user  codes  are  orthogonal,  we  can  show  that  (see  Appendix 


^B)E<3,(Iu-Ix)‘:=(l/8Tc)(^r4N+ 


(3.37) 


^Tr  (2Re(C  (1)  )+2Re{C  (1)  }+c*  ,  (-D<  IDU 

P  U  A  A  y  A  p  A 


where 


“  Cl  Q 

^T=f  ^r(s)ds=y*  Rp  (  s)  ds 


^  f  f R 


r  (s)^r  (s)  ds 


(3.38) 


(3.39) 


As  shown,  the  variance  per  interferer  depends  on  the  aperiodic 
cross-correlation  and  autocorrelation  functions  of  the  sequences 
in  C  ^  and  the  chip  correlation  functions.  For  our  preferred 
example  chip  (the  sine  pulse) 


^?r=erTc(‘293) 
^r=£rTc{  *043) 


(3.40) 


(3.41) 


This  indicates  that  the  variance  depends  only  weakly  on  the 
various  aperiodic  correlation  functions  and  could  be  approximated 


as  follows. 


2 

(6)  (S) ,T  _T  . 2^  erN 
Eb  x  (Iy V  (.  293) 


(3.42) 


Additionally,  if  our  user  codes  were  random  cosets  of  an 


.  ' .  ■ . '  .*• .  • 


.  A  ■  . 


orthogonal  code,  we  would  find 


E=°’{cx  .x'1’1"0 
y  A 


Eca){Cx  x  (J^)}=0 
C  xy 'XA 


Hence , 


E*a)  E h(B?E  (B)  (I  ~X\)2= 
c  b  t  v  X1  2T 


We  could  then  write 


"2No  2^r  ,  „ 

-*5  - 

Pr  (e)  ~  (M-l)Q 

- 2+  - L.  (K— 1 ) 

Ner  NTcer 

For  a  sine  pulse  chip  waveform 


Pr  (e)  »  (M-l)Q 

2N  cQC 

°+  ,586  (K-l) 

N  e_  N 

_ 

r  J 

- 

(3.43) 

(3.44) 

(3.45) 


(3.46) 


(3.47) 


This  result  will  be  plotted  for  various  N  and  discussed  in 
Subsection  3.4.4,  but  first  we  will  develop  a  probability  of 
error  bound  based  on  characteristic  functions. 

3.4.3  Character i s t ic  function  method 

In  this  subsection,  we  will  use  character i s t i c  functions  to 
bound  the  multiuser  error  probability  (MEP).  Characteristic 
functions  have  been  used  to  analyze  performance  in  the  presence 


of  intersymbol  interference  and  additive  Gaussian  noise  by  Shimbo 
and  Celebiler  (1971).  Additionally,  they  have  been  used  to 
analyze  the  effect  of  intersymbol  interference  or  Rayleigh  fading 
by  Vanelli  and  Shehadeh  (1974).  However,  they  also  have  been 
used  in  the  analysis  of  CDMA  systems  by  Geraniotis  and  Pursley 
(1982)  and  Geraniotis  (1983).  In  Geraniotis  and  Pursley  (1982), 
binary  and  quaternary  DS/SSMA  communications  in  the  presence  of 
additive  white  Gaussian  noise  are  considered.  The  characteristic 
function  method  is  used  to  give  probability  of  error  expressions 
when  specific  signature  sequences  (m-sequences  and  Gold 
sequences)  are  used.  In  Geraniotis  (1983),  direct  sequence  and 
frequency  hopped  SSMA  communications  are  studied.  In  this  work, 
Geraniotis  limits  his  consideration  to  random  sequences,  but  he 
includes  the  effect  of  fading  channels.  In  what  follows,  we  will 
draw  on  both  of  these  results. 

In  this  analysis,  we  will  consider  random  coset  orthogonal 
codes.  The  probability  of  error  bound  given  in  Section  3.3.4  may 
be  rewritten 

M 

Pr  (e)  <  (1/M)  I  I  Pr  (  (2/Ner)  (I.  -I  .  +  4  -vtf,  )#>1|>  sent)  (3.48) 
A=1  u^A  1  A,u  A,A  U  a  . 

M 

= ( 1/M)  I  I  (l-F(l))  (3.49) 

A=1 

where  F  is  the  cumulative  distribution  function  of  the  normalized 
noise  plus  multiple  access  interference.  If  0  (oj)  denotes  the 
characteristic  function  of  the  normalized  noise  plus 


interference,  then 


oc 


F(1)  =  (1/2tt) 


/ 


( 1/ j  oj  )  exp  (  jw)  p  (  oj  )  &u+h 


(3.50) 


Since  the  atmospheric  noise  and  multiple  access  interference 
terms  are  all  independent,  we  may  write 


K  (6) 

4>  (uj)  =<*)T  (2oj/Ne_)  7T.  <f>  (2u/Ne  ) 

■*-  ^  n _ ^  * 


(3.51) 


where  (co)  is  the  characteristic  function  of  ( I  a  —  1 A  X^  an<^ 

<f>^  (co)  is  the  characteristic  function  of  I.[0)  _  i  (0)  .  if  Equation 

H*  A 

(3.51)  is  rearranged  as  follows 


K 

4>  (^)  =4>  ( 2u)/Ner )  —  4>  ( 2u/Ne r )  (1-  7T 
n  In  1  g  =2 


<t> {B) (2w/Ner) ) 


(3.52) 


and  substituted  into  Equation  (3.49),  the  following  result  is 
achieved 

Pr  (e)  <_(M-1)  Q  (  (es/NQ)  2)  (3.53) 

oo 

M  r 

+  ( l/M)  I  Z  (1/2tt)  /  (l/ju))exp(jw)4»I(2w/Ner) 

X=1  _ro 

K 

.  (1-  TT  <$> (6)  (2oj/Ne  )  )  du> 

8  =  2 


The  first  term  on  the  right  of  this  inequality  results  from  noise 


alone  while  the  second  term  is  the  increase  in  the  probability  of 
error  due  to  the  multiple  users. 


Since  I A ^  -  I  A  ^  has  a  Gaussian  distribution,  its 
character istic  function  is  given  by 


pj  (to)  =exp  (-to  e  pNNQ/4 ) 


(3.54) 


Unfortunately,  the  characteristic  function  (Ct>)  is  not  so 

straightforward,  because,  after  all 


(6) 


)  } 


(3.55) 


where  the  and  1^0)  are  given  by  Equation  (3.25).  However  in 


Appendix  B,  we  derive 

T 

„(£) 


(to)  =  ( 1/T 


>/ 


N— 1 
7T 

n=0 


to&(x) 


J 

o  z 


/T/l-cos (2mn/N)  )  x 


(3.56) 


M-l  toR(i)  _ _ 

7T  J  (  ~~Z  / 2  /l-cos  (  27in/N )  )  dx 

n=0  °  ^ 


where  J 


Q  is  the  Bessel  function  of  order  zero 
tt/2 

J  (x)  =  (2/tt)  f  cos(xcos0)d9 


(3.57) 


The  above  expression  allows  us  to  write 


TT  *  (6)  (  2oj/Nc  r )  = 

0=2 

rp 

•  c  N-l 


(3.58) 


0)$(T) 


'((l/T  > /  *  ao,— 

c  J  n=0  ° 


JlJl  -cos(2Tm/N)  ) 


N-l 

TT 

n=0 


ujRJjO  /2  /l-cos  ( 27m/N)  )di)K  1 

'  XT- 


o'  Me 


Since  this  result  is  independent  of  /i  or  X  ,  we  get 
Pr(e)  ^  (M-1)Q(  (es/  No)is)  +  (3.59) 

°°  K 

(M-D/tt  f  (l/ai)(sino)4>I(2u/Ncr)  (1-  TT  <t>(6)  (2m/Ner))dw 

0  6=2 

Expressions  (3.58)  and  (3.59)  have  been  used  to  generate 
probability  of  error  curves  as  a  function  of  (^^/Nq)  with  K  arid  N 
as  parameters.  These  curves  are  presented  and  discussed  in  the 
next  section. 

3.4.4  Numer  ica  1  results 

In  Section  3.4.2,  we  used  the  Gaussian  assumption  to 
approximate  the  multiuser  error  probability  (MEP)  (see  expression 
(3.46)),  and  in  Section  3.4.3,-  we  used  characteristic  functions 
to  obtain  a  bound.  In  this  section,  the  numerical  evaluation  of 
the  Gaussian  and  the  characteristic  function  expressions  (Pc  and 
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PG  respectively)  will  be  briefly  discussed.  Then  the  resultant 
curves  will  be  presented  and  discussed.  Finally,  an  alternate 
character istic  function  approach  will  be  proposed  and  evaluated. 

As  shown  by  Expression  (3.46),  the  numerical  evaluation  of  PG 
is  very  straight  forward  and  the  computational  effort  is 
independent  of  K  or  N.  On  the  other  hand,  the  computational 
effort  required  for  the  evaluation  of  Pc  (inequality  (3.59)) 
does  not  increase  with  K,  but  does  increase  linearly  with  N. 

The  expressions  Pc  and  PG  are  plotted  in  Figures  3.3 
through  3.7  with  the  following  parameters: 


Figure  No. 

K 

N 

3.3 

2,3 

17 

3.4 

2,3 

31 

3.5 

3,7,11 

127 

3.6 

3,7,11 

257 

3.7 

3,7,11 

509 

All  plots  are  for  systems  where  the  chip  waveform  is  a  sine  pulse 
and  the  plots  have  the  energy  per  bit  to  noise  density  ratio  as 
the  independent  variable. 

The  characteristic  function  curves  PG  result  from  the 
numerical  evaluation  of  Expressions  (3.58)  and  (3.59).  As  such, 
they  are  approximations  instead  of  upper  bounds,  because  of 
computational  errors.  However,  we  believe  that  they  are  very 
close  to  true  upper  bounds,  because  considerable  care  was  taken 


in  the  numerical  evaluation  of  (3.58)  and  (3.59).  The  numerical 
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error  of  the  PG  curves  is  less  than  1x10“^,  which  is  acceptable 
for  the  range  of  values  plotted  in  Figures  3.3  through  3.7. 

Moreover,  we  believe  that  the  PG  curves  are  very  close  to  the 
actual  MEP  for  probability  of  error  less  than  .001.  This 
accuracy  results,  because  the  union  bound  is  very  sharp  in  this 
region . 

As  shown  in  the  figures,  considerable  discrepancy  between  the 
Pc  and  PG  curves  is  possible.  Indeed,  a  difference  of  4  dB 
exists  at  Pr(e)=4xl0-^  for  K=3  and  N=17.  A  slightly  greater 
difference  is  shown  at  Pr(e)=lxl0“3  for  K=ll  and  N=127.  However, 
no  difference  is  discernible  for  K=3,  N=509  or  K=3,  N=257  curves. 
The  discrepancy  grows  with  increasing  K  or  £b/NQ,  because  the 
relative  magnitude  of  the  multiple  access  interference  is 
increasing.  Correspondingly,  the  disagreement  shrinks  with 
increasing  N,  because  of  the  decrease  in  the  relative  magnitude 
of  the  multiple  access  interference. 

Before  concluding  this  section,  we  turn  our  attention  to  a 
third  technique  for  evaluating  the  multiuser  probability  of 
error.  This  technique  is  the  same  as  our  characteristic 
function/random  sequence  method' except  that  we  make  the 
additional  assumption  that  the  I and  I jf*'1  are  independent  for 
all  /i  not  equal  to  X.  With  this  additional  assumption,  the 
multiuser  characteristic  function  becomes 

Tc 

4>  (B)  (m)  =  (  (1  /Tc)f  (Jo(w&(T)/2)  Jq(wR(t)/2)  )Wdi)  2  (2.60) 

0 


This  new  expression  was  used  with  (3.59)  to  create  multiuser 
probability  of  error  curves  and  these  curves  were  compared  to  our 
original  Pq  curves,  which  were  not  based  on  an  independence 
assumption.  All  differences  between  the  two  sets  of  curves  were 
always  less  than  the  numerical  error  of  our  algorithm 
(Pr ( e ) <10-8)  .  We  observe  that  the  1^°)  and  are  not 

independent,  but  an  independence  assumption  does  not  produce 
meaningful  probability  of  error  discrepancies.  This  observation 
is  consistent  with  the  uncorrelatedness  of  l,\a^  and  r  (a )  for 
random  coset  codes  and  the  nearly  Gaussian  nature  of  the  multi¬ 
user  interference  pdf. 

3.4.5  Alternative  signalling  schemes 

In  the  last  few  subsections,  we  have  used  the  Gaussian 
approximation  and  characteristic  function  approach  to  provide  MEP 
expressions  for  our  CDMA  system,  which  uses  r-phase  modulation 
and  orthogonal  signal  sets.  These  techniques  can  be  used  to 
obtain  MEP  results  for  a  wide  variety  of  signalling  schemes,  and 
in  the  next  few  paragraphs  we  will  summarize  results  for  the 
following  systems: 

•  biphase  modulation/orthogonal  signal  sets 

•  biphase  modulation/antipodal  signal  sets 

•  r-phase  modulation/antipodal  signal  sets. 

Throughout  the  section,  random  coding  analysis  is  used. 

A  CDMA  system,  which  is  identical  to  ours,  except  that  it 
uses  biphase  modulation  instead  of  r-phase  modulation,  can  be 
analyzed.  For  biphase  random  coset  codes,  the  second  moment  of 
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Since  this  is  identical  to  the  second  moment  for  our  r-phase 
codes,  both  biphase  and  r-phase  modulation  result  in  the 
approximate  MEP  of  Expression  (3.46).  Consequently,  no  large 
improvement  in  the  average  multiple  access  performance  is 
realized  through  r-phase  modulation. 

Additionally,  the  MEP  for  the  biphase/orthogonal  system  can 
be  bounded  with  Expression  (3.59).  If  we  assume  that  the  1^°^ 
and  are  independent  and  that  the  chip  waveform  is  a  slow 

function  of  time  when  compared  to  the  carrier,  then 

Tc  71/2 

<J>(6)  (io) -(  C2/Tcit) y  f  (  (cos  (cu .  5^  (t  )  cos© )  )  (3.62) 

0  0 

•  '(cos  (uj.  5R(t)  cos 6)  )  )  Nd0dx)  2 

As  discussed  in  the  introductory  chapter,  a  CDMA  system  which 
uses  biphase  modulation  and  antipodal  signal  sets  has  been 
analyzed  by  Pursley  (1977),  Pursley  (1980),  and  Geraniotis 
(1982).  Under  the  Gaussian  approximation,  the  MEP  for  this 
system  is 


Pr(e)S  Q 


(3.63) 


From  Geraniotis  (1982),  the  MEP  for  this  system  that  results  from 
the  characteristic  function  method  is 


For  a  system  where  the  r-phase  modulation  is  combined  with 
antipodal  signalling,  the  Expressions  (3.63)  and  (3.64)  still 
aPPfy»  but  Expression  (3.65)  must  be  replaced  with 

Tc 

4>  (B)  (w)  =  (l/T  )f  (J  (oj$(t)/2)J  UR(T)/2))Ndr  (3.66) 

C  J  o  o 

0 


Equations  (3.58)  through  (3.66)  could  be  used  to  compare 
systems  employing  orthogonal  signalling  sets  to  antipodal 
signalling  sets.  Indeed,  the  probability  of  error  would  be  less 
for  the  orthogonal  system  for  a  given  K,  £b/No  an<^  rate. 
However,  when  rate  is  at  a  premium,  the  antipodal  system  should 
employ  an  "outer"  error  correcting  code.  This  improvement  muchly 
improves  the  performance  of  the  antipodal  system,  but  a 
comparison  incorporating  this  embellishment  is  beyond  our  current 


scope. 


the  orthogonal  system  could  also 


Additionally, 


incorporate  an  outer"  code,  where  the  outer  code  could  be  any 
variety  of  M  symbol  convolutional  or  block  code.  Once  again, 
this  improvement  is  beyond  our  current  scope.  In  this  section, 
we  have  used  character istic  functions  to  derive  the  probability 
of  error  for  CDMA  systems  using  random  sequence  sets.  In 
Appendix  C,  we  consider  the  use  of  characteristic  functions  for 
estimating  the  probability  of  error  for  deterministic  sequence 


CHAPTER  4 


SEQUENCE  DESIGN 


4.1  Introduction 

In  this  section,  we  present  two  sequence  sets,  which  achieve 
the  average  performance  of  Chapter  3  (one  is  orthogonal  and  the 
other  is  nearly  orthogonal)  and  have  bounded  worst  case  mutual 
interference.  To  do  this,  a  worst  case  analysis  of  the  multiple 
access  interference  is  presented,  and  then  a  sequence  design 
criterion  is  developed  with  the  help  of  a  result  from  number 
theory.  Next,  we  will  detour  slightly  to  consider  the  LF  multi- 
path  interference  problem  and  the  signal  acquisition  problem 
because  these  phenomena  may  impose  their  own  sequence  design 
criteria.  Finally,  we  present  some  designs  based  on  additive  and 
multiplicative  characters. 

4.2  Worst  Case  Analys i s 

The  worst  case  strategy  considers  the  interference  to  be  a 
function  of  several  deterministic  variables.  The  resulting 
sequence  design  strategy  minimizes  the  maximum  value  the 
magnitude  of  the  interference  can  take.  Accordingly,  we  seek  an 
expression  for 

l T (a, 3)  | 

max  HA  |  (4.1.) 

where  the  maximum  is  with  respect  ot  user  numbers,  data  sent,  and 
interfering  user  delays.  To  make  the  analysis  simpler  we 
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introduce  some  functions.  First  of  all,  the  correlation  between 
the  signature  sequence  x  and  the  overlap  of  the  sequences  u  and  v 
is  (same  as  Equation  3.16) 

N— 1  *  £-1  _  * 

HX,U,V{£)=  1.  xn-£un  +  Zn  Xn+N-£Vn  (4.2) 

n=£  n=0 


This  function  may  be  related  to  the  aperiodic  cross-correlation 
functions  described  in  Sarwate  and  Pursley  (1980)  as  follows 


(4.3) 


(4.4) 


where  0XfU(y)  is  the  periodic  cross-correlation  fun  lion  which  is 
also  extensively  discussed  in  Sarwate  and  Pursley  (1980).  With 
the  help  of  Definition  (4.2),  we  may  define  the  complex  multiple 
access  interference 


I ,  =^r  (t  )  exp  jut 


X  u  ,  V 


+  Rr(T)expjw  ' 


H  ( >.  + 1 ) 

x, ,u,v 

A 


(4.5) 


Contrast  this  expression  to  the  one  for  the  real  multiple  access 
interference  given  by  Equation  (3.25)  and  note  that 


I  =^Re ' I  ■  } 


(4.6) 


An  upper  bound  on  the  maximum  magnitude  of  the  real  multiple 
access  interference  may  now  be  developed  as  follows  (see  Pursley 
(1982)  ) 


max 

b,  i 

!  I j  i  <  max  h j i A ! 

b ,  t 

(4.7) 

<  max  h ( ! ft- ( x )  ]  exp j  ^  x  ; 
b ,  x 

!  H  ( i )  i .  + 

1  x ,  ,  u ,  V 

| Rr (i) II expjwcT 1 

I H  (i+1) 1 ) 

1  x> ,u,v'  1 

=max  h  (^-p  (x)  |  H  ( £ )  | 

b ,  x  1  XA'U'V 

+  Rr ( t )  ! H 

r  ■  xA,u,v 

(i  +  1)  |  ) 

observe 

that 

0  <  ( x )  <_  tr 

i  A 

(4.8) 

0  <  Rr (t)  1 


(4.9) 


0  <  1L  ( T  )  +R-,  (x)  < 

—  i  i  —  1 


(4.10) 


These  inequalities  indicate  that  the  possible  values  of 

/■v 

(Rp(T)  ,  Rp(T) )  form  a  convex  set.  Furthermore,  the  right  hand 
side  of  Inequality  4.7  is  a  linear  function  on  that  convex  set. 
C'  asequen  1 1 y ,  the  maximum  magnitude  of  that  linear  function 
occurs  at  one  of  the  corner  points  of  the  convex  set.  This  means 


max 
b ,  x 


I  1  < 
>  '  _ 


max 


u ,  V: 


max 

:C  -3) it {0, 


i,  - 
'2  -  ■ 


H 


x 


U  ,  V 


U) 


(4.11) 


N— 1  } 


For  details,  see  Pursley  (1982).  The  worst  case  interference 


given  by  Equation  (4.11)  could  be  minimized  by  minimizing  the 
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I 


1 


maximum  member  of  the  set 


{  lHx,u,v(£)  I  ;xec(a)  ;u'V£c(B) '-“^£{1,2  ...  K}; 


(4.12) 


a^8 ; { 0 , 1  ...  N— 1 } } 


A  brute  force  search  for  the  maximum  member  of  the  set  given  by 
(4.12)  would  require  the  evaluation  of  order  cross 


correlation  functions,  consequently,  we  consider  Equation  (4.3) 


and  the  Bound  (4.11)  to  derive 


b(S),T(6) 


(4.13) 


max  max  ( 

u,v  C (8)  £e{0, . . . N— 1 } 


<|Cx,u(t>  l  +  ICv,x«-f)  |) 


This  bound  is,  in  turn,  developed  through  the  use  of  a  partial 


sum  theorem,  which  we  discuss  in  the  next  section. 


4.3  Partial  Sum  Theorem 


In  this  section,  we  develop  a  partial  sum  theorem,  so  that 


we  can  bound  the  right  hand  side  of  the  inequality  (4.13).  This 


bound  will  allow  us  to  consider  "full"  sums  instead  of  the 


awkward  partial  sums  CXfU(<)  or  CVjX  (N-i. ) .  We  begin  by  defining 


a  window  sequence 


w=  { w  } 


(4.14) 


1  ms {0 , 1 . . .N-l } 


0  otherwise 


"  -  ",  "4  JV  - 


w  '  w  « 


and  by  considering  the  sequences  x  and  u  to  be  infinite  length 
sequences  with  periodicity  N.  Next,  the  associated  ambiguity 
functions  are  defined  as  follows 


A  ( i , c) 
w ,  w 


=  I  w  w  exp  ( j  2TTcm/N) 
„  m  m+i. 
m=0 

N-Ji-1 

=  I  exp ( j 2ncm/N) 

m=0 


(4.15) 


4X/U(H.c)  =  I_o  Xnun+Jt  exp(j2ircn/N) 


(4.16) 


Observe  that 


A  (l,c)  A  (£,c)  = 

Q  x ,  u  v  w ,  w 

1  *  N-l  N-l 

,  xnun+l  *  wmwm+  a  *  exp ( 32*0 (n-m> /N)  = 
j  m=0  c—  U 

N-l 

T  x  w  u  ,,w  ,0  =  N  C  (l)  = 
n  n  n+S.  n+x.  x,u 

n-u 

Vu(t'0)  +  Vu**'11  Vw(1'°>  + 


(4.17) 


+  4X,U(1'N-l>  Vw'1'*'11 


I  <l)  (N-l)  +  X  4  <1,0  AW(Un.,o 

A  ,  U  c=: 


Consequently , 


NC  (l) 
X  ,  u 


(N-fc)  |  9  (  2.)  | 

—  1  v  _  n  1 


(4.18) 


+  max  [A  n(ifC)  |  I  1 A  ..U,c) 
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where  [N/6 ]  denotes  the  greatest  integer  less  than  N/6.  We  have 
computed  the  left  hand  side  of  Equation  (4.21)  for  N  between  20 
and  1000  and  find  that  in  this  range,  we  have 

N-l  |  sin2cTr/N  | 

( 1/N)  max  I  -  <_  .  44  751nN  N  _>  20  (4.22) 

2  c=l  |  sinnc/N  j 

The  coefficient  of  InN  decreases  from  .4475  as  N  increases  from 
N  =  2  0  to  N  =  1 0  0  0 .  Consequently,  we  use  the  coefficient  .4475  in 
what  follows  and  we  feel  that  this  bound  is  valid  for  all  N  .cp 
12  greater  than  20.  Combining  these  results  with  the  Inequality 
(4.18)  yields 

max  !cxn^)limax  {|0V,,U)|  (N-JU/N)  +  (4.23) 

x,u,2  x'u  x , u, 2  X'U 

•  .4475  InN  max  { | A  ( 2 , c ) I } 

x,u,2,c^0  x'u 

Additionally,  we  may  now  rewrite  Inequality  (4.13)  as 

max  |H  (2) [  <  max  je  (2)|  +  (4.24) 

X,U, v, 2  '  '  x , u , 2  ,u 


.395  InN  max  |A  (2,c) I 

x,u,2,c^0  x,u 

This  approach,  which  is  similar  to  an  approach  described  by 
Lerner  (1961),  allows  us  to  consider  the  function  A  (  (  ,c) 

A  /  LI 

instead  of  the  more  difficult  partial  sum  Cx  u(*.  ).  Before 
turning  our  attention  to  the  presentation  of  sequence  set 
designs,  we  describe  the  effect  that  the  LF  multipath  problem  and 
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the  signal  acquisition  problem  have  on  our  sequence  design 
strategy. 


4.4  Low 


Multipath  and  Sequence  Design 


For  most  of  the  LF  band,  multipath  is  specular  in  nature  and 
depends  on  the  transmitter  to  receiver  distance.  It  is  a  very 
important  part  of  any  LF  channel  model.  For  example,  at  500  km, 
the  multipath  signal  (skywave)  can  exceed  the  groundwave  signal 
strength  by  tens  of  decibels  at  100  kHz.  On  the  other  hand,  at 
100  km,  the  groundwave  is  generally  at  least  10  dB  stronger  than 
the  skywave.  For  long  ranges,  skywave  is  problematic,  because 
the  delayed  signal  may  resemble  a  signal  other  than  the  one 
currently  being  communicated.  Using  an  analysis  similar  to  the 
worst  case  multiple  access  analysis,  it  can  be  shown  that  the 
worst  case  skywave  delay  is  equal  to  an  integer  number  of  chip 
times  ( )  .  In  such  a  situation,  the  interference  introduced 


by  skywave  in  the  receiver's  correlator  is 

h£-Re{H  (\  T  ) } 

i  x,u,v  ^s  c 


(4.25) 


where  x  =  x^a^  and  u  and  v  are  the  previous  and  the  current 
sequences,  respectively,  being  sent  by  the  transmitter. 
Since  the  skywave  delay  will  always  be  small,  compared  to  the 
duration  of  any  sequence  transmission,  the  magnitude  of  the 
skywave  interference  may  be  conveniently  bounded  as  follows 

’>M1!x,u.v<i's’  •  r  2  's' 


Skywaves  of  appreciable  strength  compared  to  groundwave  are  not 


delayed  by  more  than  200  /isec,  which  for  our  model  means  that  K  s 
will  not  exceed  3.  Hence,  the  maximum  value  of  the  following  set 
is  a  reasonable  skywave  sequence  design  criterion 

{  I  6x,v(^s)  I  ;x'vec(a)  ;  { 1 , 2 , 3 ) ]  (4.27) 

4.5  S ignal  Acquisition  and  Sequence  Design 

Signal  acquisition  is  the  process  of  making  the  correlators 
in  the  a*-*1  receiver  time  synchronous  with  the  reception  of 
the  a*-*1  signal.  For  our  phase  coherent  receivers,  this  process 
may  be  described  as  a  3  part  procedure,  where  the  parts  are 
carrier,  chip,  and  code  word  synchronization.  Carrier 
synchronization  is  the  process  of  bringing  the  receiver's  carrier 
reference  into  phase  with  the  received  carrier.  Correspondingly, 
chip  synchronization  is  the  procedure  the  receiver  uses  to  align 
its  chip  reference  time  with  the  received  chip  waveforms.  Many 
carrier  and  chip  synchronization  schemes  are  described  in 
Stiffler  (1971).  Once  carrier  and  chip  synchronization  have  been 
achieved,  the  problem  of  code  word  or  sequence  synchronization 
remains.  Code  word  synchronization  is  the  problem  of  identifying 
the  first  element  of  the  received  sequences,  or  equivalently, 
determining  whether  the  currently  observed  string  of  elements  is 
a  single  complete  sequence  of  an  overlap  of  two  sequences.  If  no 
external  synchronization  aid  is  available,  this  process  may  be 
enhanced  by  using  comma  codes  or  comma  free  codes.  In  this 
paper,  comma  codes  will  be  considered. 


In  our  case,  comma  coding  will  mean  the  transmission  of  a 
known  member  of  the  sequence  set  (the  comma)  at  known  intervals. 
Such  a  comma  is  called  singular  if  the  comma  is  guaranteed  to  be 
different  from  any  code  word/code  word  overlap.  However,  for 
noisy  channels,  comma  coding  schemes  do  not  usually  involve  this 
constraint.  Instead,  the  receiver  sums  the  output  of  the 
correlator  matched  to  the  comma  over  several  comma  repetition 
intervals,  and  this  sum  becomes  the  decision  statistic  for  the 
code  word  synchronization  algorithm.  If  the  comma  is  being 
observed,  the  mean  of  the  sum  increases  linearly  with  each 
repetition  interval.  If  a  code  word/code  word  overlap  is  being 
observed,  the  ratio  of  the  sum's  standard  deviation  to  the 
synchronized  mean  can  be  made  arbitrarily  small  with  increasing 
observations.  If  a  comma/code  word  overlap  is  being  observed, 
the  test  also  becomes  increasingly  reliable,  but  only  if  the 
comma  is  easily  distinguished  from  shifted  versions  of  itself. 
This  requirement  may  be  met  through  the  minimization  of  the 
following  sequence  design  criterion 


max  min  {  1  cx  ^  ^  )  !  ;xeC(ct)  ;U(1,2.  .  .  N  - 1 J } 


(4.28) 


This  criterion  guarantees  that,  at  least,  one  member  of  each 
sequence  set  is  a  suitable  comma.  In  our  analysis,  bounds  for 
the  worst  case  Cx  x  (.)  in  every  code  will  be  found  by  applying 
the  partial  sum  theorem  in  a  manner  similar  to  the  multiple 
access  analysis.  We  may,  finally,  turn  our  attention  to  the 
description  of  specific  sequence  sets. 
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4.6  Sequence  Designs  Based  on  Additive  Characters 

In  this  section,  sequence  designs  based  on  additive 
characters  are  considered.  We  will  define  additive  characters, 
describe  some  of  their  fundamental  properties,  and  state  a  useful 
additive  character  sum  theorem.  Next,  we  describe  some  additive 
character  designs,  which  have  been  presented  in  the  literature. 
Finally,  we  propose  a  sequence  design  for  our  LF/CDMA  model  and 
analyze  its  multiple  access,  signal  acquisition  and  skywave 
performance . 

An  additive  character  ^  is  a  map  from  the  additive  group  of 
GF  (p)  co  the  complex  numbers  of  unity  magnitude,  such  that  for 
all  x  and  y  in  GF  (p) 

<1*  (x+y)  (x)  ip  (y)  (4.29) 

This  definition  implies  that  ^{0)=1  and  that  every  additive 
character  on  GF(p)  is  of  the  form 

ip  (x)  =exp  j  2Trax/p  (4.30) 

for  some  integer  a.  This  expression  is  revealing,  because  it 
shows  that  sequence  designs  based  on  additive  characters  will 
employ  p  phases.  Thus  p  should  not  be  too  large,  because  of 
implementation  considerations.  However,  the  following 
exponential  sum  theorem  (Schmidt  (1976))  allows  us  to  prove  some 
remarkable  correlation  results  for  these  sequences.  Let 

g(x)=a  xn  +  a  ,xn  ■*■+... a.  (4.31) 

n  n-1  0 


be  a  polynomial  with  integer  coefficients  and  0<n<p  and  p  lan. 
Then 

|  X)  expj2rrg  (x) /p  j  £  (n-1)  /p~  {4  32) 

x£GF(p) 


Additive  character  sequences  of  period  p  (N=p)  are 
considered  in  Frank  and  Zadoff  (1973)  and  Chu  (1972).  If  we 
limit  consideration  to  their  sequences  of  odd  prime  period,  the 
kt*ri  element  of  one  of  their  sequences  is  given  by 

vk=  ^(ak2)  ae{l  ,2. .  .p-1} 


These  sequences  have  the  ideal  periodic  auto-correlation  function 


9  (5  ) 

v, v  K 


N  51  =0  mod  N 
0  otherwise 


(4.34) 


Also,  the  sequence  definition  may  be  generalized  by  the  addition 
of  a  linear  phase  shift  term  without  changing  the  auto¬ 
correlation  behavior  (Chu  (1972)). 

2  ^ 

vk=i|;1  (ak  +bk)  be  { 0 , 1 .  .  .p-1 }  (4.35) 


Later  on,  sets  of  sequences  were  defined  and  analyzed  by 
Sarwate  (1979)  and  Alltop  (1980).  In  their  case,  the  element 

of  the  a*"*"1  sequence  is  given  by 


v 


ae { 1 , 2 . . . o-l 


(4.36) 


This  set  contains  p-1  sequences  of  period  p  (N  =  p)  and  the 
periodic  correlation  function  obeys 


M  v?v' 

'  N  v=v  ‘,1=0 

0  v=v  1  j£0 


(4.37) 


Additionally,  sets  of  cubic  phase  sequences  were  considered  by 
Alltop  (1980),  who  defined 

v  =p  (k3+ak)  ae { 1 , 2 . . . p-1 }  (4.38) 

3.K  1 

These  sequences  are  also  periodic  p,  but  their  periodic 
correlation  functions  have  slightly  different  character ist ics 


rt) 


N  v=v 1  , 5  bO 
•  0  Vj*V ',5=0 

/iT  otherwise 


(4.39) 


A  sequence  set  for  use  with  our  LF/CDMA  model  is  now 
proposed.  Let  the  Afc^  sequence  (data  =  A)  in  the  code  C  ^  be 
defined  by 


p-1 

} 

k  =  0 


ac{l,2...p-l;  >  e  r  0 , 1  ...  p-1}  (4.40) 


where 


k) 


(4.41) 


This  definition  orovides  p-1  codes  (K<p-1)  with  p  sequences  (M=p) 


r 


in  each  code  and  every  sequence  has  length  or  period  p  (N=p). 
Additionally,  each  user's  code  is  a  complete  orthogonal  set, 


because  for  all  v  and  v' 


6  ,  (0)  = 
v  v 


0  v^v’ 


N  v-v ' 


(4.42) 


For  the  analysis  of  our  design  with  respect  to  multiple 
access  performance  and  comma  coding  perfomance,  we  consider  the 
following  function  as  recommended  in  Sections  4.3  and  4.4. 


,  (S,< 


YL  Vk[vk  +  \]  expj2irck/N  j 
k=0  K  K+x 


(4.43) 


Y  expj2m  (k3  (a- a  ' )  +k2  (-a  '  3^  )  +k  ( A-3a  '  ,l2-  '  +c)  )  /o’ 
ksC-F  (p) 


For  the  multiple  access  case  {a^Ct') ,  we  apply  the  exponential  sum 
theorem  and  find 


' Av,v'  (1,c)  Y2' 


(4.44) 


We  continue  the  analysis  by  applying  the  partial  sum  theorem  as 
f o 1  lows 


H  (  J)  |  <  2/d  +  1.79  0  InM  ^ 

X  ,  U  ,  V  '  ’  ' 


(4.45) 


This  bound  may,  in  turn,  be  used  to  approximately  bound  t h • 
maximum  interference  due  to  a  single  competing  user  as  follows 


u ,  v. 


U  (  O 

X ,  U  ,  V 


M  (  1 + . 39  5  InN ) 


(4.4b) 


% 


This  result  does  not  guarantee  that  the  worst  case  interference 
is  less  than  the  signal  strength  unless  sequences  of  length 
exceeding  110  are  used.  Consequently,  when  considering  multiple 
competing  users,  Equation  (4.46)  cannot  be  used  to  bound  the 
user's  probability  of  error  in  any  reasonable  way.  However,  it 
does  show  that  the  worst  case  multiple  access  interference  to 
signal  strength  ratio  can  be  made  arbitrarily  small  with 
increasing  N  and  does  provide  a  way  for  comparing  different  code 
set  designs. 

To  consider  the  signal  acquisition  performance  of  our 
design,  we  apply  the  exponential  sum  theorem  to  Equation  (4.43) 
with  a=a'  and  A=A'  and  find 

(4-47) 

Applying  the  partial  sum  theorem  yields 

i  -  4475  /N  lnN  (4.48) 

This  result  is  somewhat  disappointing,  because  it  shows  that  N 
must  equal  110  before  the  worst  case  auto-correlation  sidelobe  to 
sequence  length  ratio  equals  0.2.  However,  it  does  show  that 
with  increasing  N,  every  sequence  in  every  code  can  serve  as  .a 
comma.  The  bound  of  (4.48)  suggests  that,  at  least,  one  very 
good  comma  can  be  found  in  each  code,  if  some  additional  search 
scheme  is  designed. 

As  mentioned  in  Section  4.4,  the  s  k  y  wave  rejection 
performance  of  any  single  set  may  be  bounded  through  e’-a 1  uat  i  on 


of  the  worst  case  magnitude  of  the  periodic  cross-correlation 
function 


ev,V<&>Hl=o  vk[v'+)J*| 


=  |  I  expj27T(k2(-a'35l}+k(X-3a^  2-X')  ) /d| 
keGF(p) 


Applying  the  exponential  sum  theorem  for  8^0,  we  find 


(4.49) 


max  |0vv.(l)|£^“ 
8 , v,v* 


(4.50) 


This  inequality  along  with  (4.26)  shows  that  the  worst  case 


skywave  interference  is  bounded  by 


v  Jw111  -  ^4i/3+6) 

X  /  u  /  V  f.\ 


(4.51) 


If  the  received  skywave  power  equals  the  received  groundwave 
power,  then  the  demodulated  skywave  to  groundwave  voltage  ratio 
will  not  exceed  0.2,  if  the  sequence  length  exceeds  25. 

Our  additive  character  sequence  design  ((^-design)  has  now 
been  evaluated  with  respect  to  worst  case  multiple  access,  comma 
coding  and  skywave  performance.  It  can  be  compared  to  other 
designs  through  these  bounds,  but  the  greatest  liability  of  the 
design  may  be  that  it  requires  N  phases.  Since  this  requirement 
will  make  it  difficult  to  realize  designs  with  long  sequences,  we 
turn  our  attention  to  designs  based  on  multiplicative  characters. 


This  result  does  not  guarantee  that  the  worst  case  interference 
is  less  than  the  signal  strength  unless  sequences  of  length 
exceeding  110  are  used.  Consequently,  when  considering  multiple 
competing  users,  Equation  (4.46)  cannot  be  used  to  bound  the 
user's  probability  of  error  in  any  reasonable  way.  However,  it 
does  show  that  the  worst  case  multiple  access  interference  to 
signal  strength  ratio  can  be  made  arbitrarily  small  with 
increasing  N  and  does  provide  a  way  for  comparing  different  code 
set  designs. 

To  consider  the  signal  acquisition  performance  of  our 
design,  we  apply  the  exponential  sum  theorem  to  Equation  (4.43) 
with  a=a'  and  X=A'  and  find 


Av,v{5'c) 


(4.47) 


Applying  the  partial  sum  theorem  yields 


C  (2.)  I  <  .  4475  /N  InN 
v,v  1  — 


(4.48) 


This  result  is  somewhat  disappointing,  because  it  shows  that  N 
must  equal  110  before  the  worst  case  auto-correlation  sidelobe  to 
sequence  length  ratio  equals  0.2.  However,  it  does  show  that 
with  increasing  N,  every  sequence  in  every  code  can  serve  as  a 
comma.  The  bound  of  (4.48)  suggests  that,  at  least,  one  very 
good  comma  can  be  found  in  each  code,  if  some  additional  search 
scheme  is  designed. 

As  mentioned  in  Section  4.4,  the  skywave  rejection 
performance  of  any  single  set  may  be  bounded  through  evaluation 


4.7  Sequence  Design  Based  on  Multiplicative  Characters 

As  mentioned  at  the  end  of  the  last  section,  multiplicative 
character  sequences  are  now  considered,  because  they  will  require 
fewer  than  N  phases.  In  this  section,  we  begin  by  defining 
multiplicative  characters,  discussing  their  properties,  and 
stating  two  character  sum  theories  that  will  allow  the 
calculation  of  correlation  bounds.  Then,  we  discuss  sequence 
sets,  which  have  appeared  in  the  literature,  and  finally,  we 
propose  a  multiplicative  character  sequence  design  for  our  LF 
CDMA  model. 

Let  GF(p)  denote  the  finite  field  of  order  p  and  GF*  (p) 
denote  the  multiplicative  group  of  GF(p),  consisting  of  the  non¬ 
zero  elements  of  GF(p).  Additionally,  let  g  denote  a  generator 
for  GF* (p) ,  such  that 


I  (x) 

x=g 


(4.52) 


for  all  x  in  GF  (p)  and  where  I(x)  is  the  index  of  x.  The 
multiplicative  character  on  GF* (p)  is  a  map  X  from  GF*(p)  to 
the  complex  numbers  with  unity  magnitude  such  that 


X (xy) =X (x) X (y) 


(4.53) 


for  all  x  and  y  in  GF  (p).  This  definition  implies  that  X  ( 1 )  =  1 
and  that  every  multiplicative  character  on  GF* (p)  is  of  the  form 


X  (x)  =expj27ral  (x)  /  (p-1) 


(4.54) 


Every  character  will  have 


p-1 

*a  =x0 


where  XQ  is  the  principal  character,  which  means  XQ(x)  =  l  for 
all  x.  The  order  of  Xa  is  the  smallest  positive  integer  d  such 
that  Xad=  X0.  The  order  of  Xa  equals  (p-1)  /gcd  (a, p-1) . 

The  order  is  a  particularly  important  parameter  in  sequence 
design,  because  it  is  equal  to  the  number  of  phases  used  by  the 
sequence.  The  definition  of  multiplicative  characters  may  be 
extended  as  follows  to  map  all  of  GF(p)  to  the  complex  plane 


(4.55) 


x(Q)  = 


(4.56) 


X^X0 


Additionally,  the  following  identities  hold 

X*(x)  =  [Xa(x)]P“2  (4.57) 

Xa,  (x)  =  [xa(x)  3  a 


where  a'=*aa  mod  p-1.  Two  results  on  character  sums,  which  will 
help  us  find  correlation  results,  are  now  presented.  The  first 
of  these  is  due  to  Sidel'nikov  (1969). 


I  X  t (x+c )  (x+c')  ] 

xeGF(p)  a 


(4.58) 


where  0<i<d,  c  ^  c'  and  d  is  the  order  of  Xa-  The  second  result 


is  presented  in  Schmidt  (1976)  and  bounds  the  magnitude  of  a 
multiplicative  character  sum.  Let  X  be  a  multiplicative 
character  of  order  d>l,  and  let  f(x)  be  an  element  of  GF(p)[x] 
(GF(p)[x]  is  the  ring  of  polynomials  with  coefficients  from 
GF  (p) .)  with  m  distinct  zeroes.  Additionally,  f(x)  is  not  a  dth 
power,  which  means  it  cannot  be  expressed  as  f(x)=c(  k  (x))d  where 
c  is  in  GF  (p)  and  X(x)  is  in  GF(p)  [x]  ,  then, 

I  I  X(f(*))|  1  (m-l)/p  (4.59) 

xeGF(p) 

We  now  describe  some  multiplicative  character  sequence  designs, 
which  have  been  presented  in  the  literature. 

Length  p-1  sequences  are  described  in  Lerner  (1961)  and 
similar  sequences  of  any  length  are  defined  by  Scholtz  and  Welch 
(1978).  We  will  describe  the  length  p-1  sequences  only,  because 
the  value  of  the  extension  for  our  LF/CDMA  model  is  not  clear.  A 
sequence  of  length  p  from  a  set  of  p-2  sequences  may  be  defined 
as  follows 

•  p-1 

v  ={v  ,  }  ae{l,2. . .p-1}  (4.60) 

a  aJc  k=0 

where  valc  =  Xa(X)  and  va0=0.  The  periodic  auto-correlation 
function  for  any  of  these  sequences  is 

p-1  SisO  mod  p  (4.61) 

-1  1^0  mod  p 


The  magnitude  of  the  periodic  cross-correlation  function  is  given 
by 


Vv<f>i  = 


/p 

0 


?  ^0  mod  p 
j =0  mod  p 


(4.62) 


This  sequence  set  is  nearly  ideal  with  respect  to  simultaneous 
minimization  of  the  maximum  value  of  the  periodic  auto¬ 
correlation  and  cross-correlation  functions  (Sarwate  (1979)). 

A  single  sequence  design,  which  has  a  nearly  ideal  periodic 
auto-correlation  function  is  given  by  Sidelnikov  (1969). 
Limiting  our  interest  to  GF(p),  Sidelnikov's  sequences  are 
defined  as  follows 

(4.63, 

x(gk+l)  gk+l^iO  mod  p-1 

vk=  , 

1  g*+l=0  mod  p-1 

These  sequences  have  period  p-1  (N=p-1)  and  vk  is  defined  to  be  1 
k 

when  g*+l  is  congruent  to  zero  to  provide  uniform  envelope 
signals.  Sidelnikov  proves  that 

|  S  (<J  |  < 4  )  ^0  mod  N  (4.64) 

Additionally,  if  the  order  of  the  character  is  2  (providing 
biphase  sequences) ,  then 


l9v.v(*>l=42 


t ^0  mod  N 


(4.65) 


A  sequence  set  design,  based  on  multiplicative  characters. 


(Krone  and  Sarwate  (1982))  uses  the  definition 


v  ={v 

c  c,k  k=0 


ce{0,l...p-l} 


(4.66) 


X(g^k+2gk+c)  if  g^k+2gk+c^0 


if  g^k+2gk+c=0 


This  definition  provides  p  sequences  of  period  p-1  and  it  may  be 


shown  that 


Vv(<>l 


<  3/d  +  5 


=  p-1 


J IjiQ  or  c^c' 


and  c=c ' 


(4.67) 


This  design  is  powerful,  because  the  number  of  phases  employed  by 
all  p  sequences  is  equal  to  the  order  of  the  character.  Hence,  a 
sequence  set  for  quadriphase  signalling  could  be  designed  by 
choosing  a  and  p  such  that  gcd (a ,p-l)  =  (p-1) /4. 

We  now  propose  a  sequence  design  for  our  LF  CDMA  model  based 
on  multiplicative  characters.  We  show  that  each  user's  set  is 


nearly  orthogonal  and  then  analyze  multiple  access,  signal 
acquisition  and  skywave  performance.  Let  the  Xth  code  word  of 
C^  be  defined  by 


P-2 

} 

k=0 


cte  {1 , 2  .  .  .p-2  }  Xe { 1 , 2  .  .  .  p-2  }  (4.68) 


where 

VX?k=^(gk+A)  (4*69) 

This  design  provides  p-2  sets  of  sequences  (K=p-2) ,  where  each 
set  contains  p-2  sequences  (M=p-2)  of  length  or  period  p-1 
(N  =  p-1).  We  do  not  make  an  effort  to  insure  that  every  element 
has  unity  magnitude,  because  uniform  envelope  signals  are  not 
necessarily  required  for  the  LF  channel. 

As  shown  by  the  definition,  the  user  number  (Of)  determines 
the  order  of  the  character  and,  consequently,  the  number  of 
phases  required  by  that  user.  Consequently,  if  p-2  users  (K=p-2) 
were  operating,  then  ^(p-l)  of  these  would  use  p-1  phases  where 
01  is  Euler's  totient  function.  These  J2f(p-1)  users  would  be  using 
N  phases  and  it  would  appear  that  the  design  had  failed  to 
overcome  the  main  drawback  of  additive  character  designs. 
However,  all  CDMA  systems  must  have  much  fewer  than  N  active 
users,  based  on  average  performance  analysis.  For  those  systems 
without  large  numbers  of  inactive  potential  users,  the  sequence 
length  can  be  chosen  such  that  all  users  can  be  accommodated 
without  requiring  any  to  use  a  large  number  of  phases.  For 
example,  if  we  choose  p=257  (and  consequently  N=256)  the 


user/phase  breakdown  for  31  users  could  be 


number  of  users 


number  of  phases 


In  general,  the  choice  of  the  user  numbers  (a)  must  be 
constrained  in  a  special  way  to  guarantee  good  multiple  access 
performance.  Our  design  requires  that  a  be  even  and  that 
if  a=2  a'  mod  p-1  for  any  two  codes  a  and  a',  then  ct  is 


even. 


We  may  show  that  each  individual  signal  set  C is  nearly 


orthogonal  as  follows 

ev,V(0)=I=0 

=  Z  xaMgk+M  (gk+A  ’ ) d_1  ]  -v  ,  [v '  ] 

keGF(p)  a  p  1  p 

Using  Sidelnikov's  results 


(4.70) 


ev,V«°'-1-vp-llvp-l1 


(4.71) 


Hence, 


e„  ...  (0)  I  <2 


v,v’ £C (a) ;v/v' 


(4.72) 


We  now  evaluate  the  designs  multiple  access  performance  by 
using  the  function 

N-l 

|Av  v,  (V  ,c)|=  I  z  vkv*+^expj27Tck/N|  (4.73) 

k  “  0 

ill  X„,  l  (gk+X)a(gk+X'g"^)  (gk)c] |  +1 

k£GF(p) 


where  a-a  a'  mod  p-1  and  c  =  ca'  mod  p-1.  Applying  the  sum  theorem 
for  ai 1,  it  can  be  shown  that 


max  |A  ,  ( H,c)|i2/p  +1 

ce{l,2. . .N-l}  ' 


(4.74) 


|Av,v*  U,0)|</p  +1 


Consequently, 

max  Her  |H  (i)  |  <  J5er((/N+T  +l)  +  .895  InN  (2/N+I  +1))  (4.75) 

u,v,l  1  X,U,V  “  1 

This  final  bound  on  the  maximum  multiple  access  interference  for 
our  multiplicative  character  sequence  set  (X  set)  is 
approximately  the  same  as  the  bound  for  our  ^  set. 

We  analyze  the  signal  acquisition  of  the  X  set  by 
setting  a  =1  and  A=  in  Equation  (4.73)  and  again  applying  the 
sum  theorem,  which  gives  us 

&VfV'  (  /^,c^0)|  1  2/p  +1  XfO  (4.76) 

=  P 


A..  (\'  ,  C  =  0 )  I  <  vp  +1  {  aO 


These  bounds  yield 


|C*  U)  |  <  (/N+l  +1)  +  .4475  InN  (2/N+I  +1)  (4.77) 

v  f  V  *” 

This  performance  bound  is  approximately  twice  that  obtained  for 
our  additive  character  design,  but  the  same  discussion  applies. 
Finally,  to  investigate  skywave  performance 

lev  v'^l  III  Xa,  [  (gk+A)  (gk+X'g-t)p-2]  |+1  (4.78) 

'  keGF(p)  a 

Applying  the  sum  theorem  yields 

IVv^l  (4.79) 

<N  X=X ' 

As  shown,  the  average  skywave  interference  for  our  X  design 
cannot  be  reasonably  bounded,  but  each  code  may  be  altered  as 
follows  to  provide  good  skywave  performance.  For  each  x ^  in 
C  expurge  the  set 

{x^a) ;y=Ag-S  £e{l, 2  ...  fcg}}  (4.80) 

The  new  code  will  have  worst  case  skywave  interference  bounded  by 
(1/2)  h/ff  +6),  but  it  also  has  reduced  rate. 

In  this  section,  a  design  that  can  employ  many  fewer  than  N 
phases  has  been  discussed.  In  the  final  chapter  of  this  thesis, 
we  will  summarize  these  results. 
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CHAPTER  5 

SUMMARY  AND  CONCLUSIONS 

In  this  thesis,  we  have  proposed  and  analyzed  a  code 
division  multiple  access  communication  system,  which  is 
especially  suited  for  the  low  frequency  channel. 

Our  CDMA  scheme  is  similar  to  classical  CDMA  systems  in  the 
following  respects.  It  has  K  users  sharing  a  channel  by  phase 
modulating  their  transmissions  with  signature  sequences.  These 
users  make  no  attempt  at  frequency  separation  and  the  scheme  is 
asynchronous,  which  means  that  the  users  are  not  time  coordinated 
in  any  way. 

However,  our  scheme  is  more  complicated  than  the  classical 
systems  in  the  following  respects.  Each  user  has  a  sequence  set, 
which  consists  of  M  orthogonal  sequences.  This  change  increases 
the  complexity  of  the  system,  but  now  log2M  bits  of  information 
are  transmitted  by  choosing  among  the  signature  sequences. 
Moreover,  for  a  given  bandwidth  and  communication  reliability, 
the  information  rate  is  greater  than  that  achieved  with 
traditional  antipodal  sequence  sets. 

Additionally,  our  scheme  employs  r-phase  modulation  (the 
signature  sequences  are  r-valued).  This  change  also  represents 
an  increase  in  complexity,  but  probably  provides  an  improvement 
in  the  worst  case  multiple  access  performance  of  the  system. 
Indeed,  both  of  the  sequence  designs  we  found  for  our  LF/CDMA 
communication  system  involved  r-valued  sequences. 


.V  A  .V  S' 


Our  analysis  did  not  model  each  user's  initial  carrier  phase 
as  random,  because  LF  transmitter  techniques  are  such  that  the 
initial  carrier  phase  is  fixed  and  stable  with  respect  to  the 
signal  envelope.  However,  if  the  initial  carrier  phase  was 
modelled  as  a  uniform  random  variable  on  [0,27),  our  system 
analysis  results  would  not  change.  This  invariance  exists, 
because  we  assumed  the  modulation  phase  of  each  chip  was  uniform 
on  [0,27).  The  invariance  would  also  result  for  biphase 
modulation  if  the  chip  envelope  varied  slowly  compared  to  the 
carrier . 

Finally,  our  model  includes  an  additive,  impulsive  noise 
source  to  represent  LF  atmospheric  noise.  We  used  a  truncated 
Cauchy  probability  distribution  function  to  represent  the  first 
order  pdf  of  this  atmospheric  noise  after  filtering  and  sampling. 
The  strongly  non-Gaussian  nature  of  these  noise  samples  leads  us 
to  derive  and  analyze  a  locally  optimum  Bayes  detector  (LOBD) , 
which  is  based  on  our  noise  model.  The  derived  receiver  consists 
of  a  filter  matched  to  the  carrier  modulated  by  the  chip 
waveform,  followed  by  a  sampler.  The  samples  are  then  processed 
by  a  zero  memory  nonlinearity  and  then  fed  to  a  bank  of  M  matched 
filters,  where  each  filter  corresponds  to  one  of  the  M  members  of 
the  sequence  set.  The  performance  of  this  structure  is  a  strong 
function  of  Fisher's  information  number,  which  in  turn  depends 
entirely  on  the  statistics  of  the  sampled  noise.  By  evaluating 
Fisher's  information  number  for  our  combined  Gaussian,  Cauchy  and 
multiple  access  noise,  we  found  that  every  LF/CDMA  receiver 
should  employ  a  ZNL  to  limit  the  effects  of  the  impulsive  events. 
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We  also  found  that  such  a  nonlinear  receiver  performs 
approximately  as  well  as  a  linear  receiver  operating  in  an 
environment  without  the  Cauchy  noise. 

Consequently,  we  analyzed  the  single  and  multiuser 
performance  of  a  linear  receiver  with  additive  white  Gaussian 
noise.  The  multiuser  analysis  is  approached  through  a  Gaussian 
assumption  and  a  characteristic  function  method.  The  former 
technique  is  based  on  the  assumption  that  the  multiuser 
interference  has  a  Gaussian  pdf  and  leads  to  an  approximate 
multiuser  error  probability  (MEP)  expression.  The  approximation 
is  conveniently  simple,  but  our  characteristic  function  analysis 
indicates  that  it  does  not  always  give  accurate  results.  The 
difference  between  the  results  of  the  two  approaches  increases 
with  the  magnitude  of  the  multiple  access  noise  relative  to  the 
atmospheric  noise.  The  results  of  both  analyses  are  summarized 
in  Figures  3.3  through  3.7. 

The  multiple  access  performance  of  our  LF/CDMA  model  can  be 
appreciated  by  considering  the  following  example.  Suppose,  we 
wish  to  achieve  a  MEP  of  approximately  10-5  and  we  intend  to  use 
chip  waveforms  of  100  sec  duration,  which  corresponds  to  a 
bandwidth  of  approximately  20  kHz.  Then  a  scheme  with  sequences 
of  length  17  will  accommodate  2  users  and  an  energy  per  bit  to 
noise  density  ratio  of  no  less  than  16  dB.  Each  user  will  be 
able  to  communicate  at  an  information  rate  of  2400  bits  per 
second.  On  the  other  hand,  a  scheme  with  N=257  will  accommodate 
3,7,  or  11  users  with  eb/N0  9reater  than  or  equal  to  6  dB,  8  dB, 


<*  .-  -■  S'.  S .  . 
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or  12  dB  respectively.  The  information  rate  would  be  slightly 
greater  than  300  bits  per  second. 


In  Chapter  4  of  the  thesis,  we  developed  a  sequence  set 
design  technique  for  our  LF/CDMA  system  and  used  it  to  synthesize 
2  sequence  sets. 

The  sequence  set  design  strategy  is  to  minimize  the  maximum 
magnitude  of  the  interference  that  any  one  user's  transmitter  can 
introduce  into  the  receiver  of  another  user.  Accordingly,  we 
bound  the  maximum  value  this  interference  can  take  as  a  function 
of  the  interferer's  delay  and  choice  of  sequences.  However, 
evaluation  of  this  initial  bound  requires  consideration  of  all 
the  aperiodic  cross-correlation  functions  between  the  sequences 
of  the  two  users.  Such  an  evaluation  by  computer  would  be 
prohibitively  time  consuming  and  direct  analysis  is  difficult, 
because  aperiodic  cross-correlation  functions  are  partial  sums, 
for  which  few  results  exist  in  the  literature.  Consequently,  we 
considered  a  partial  sum  theorem  (Vinogradov,  1949),  which 
allowed  us  to  bound  the  magnitude  of  partial  sums  in  terms  of  the 
magnitude  of  the  more  convenient  and  widely  studied  full  sums 
(Vinogradov  (1949),  Lerner  (1961)).  As  such  ,  we  developed  a 
multiple  access  sequence  design  criterion,  which  depends  on 
analytically  tractrable  full  sums  (ambiguity  functions). 

Even  though  multiple  access  interference  is  the  main  concern 
of  the  thesis,  skywave  (specular  multipath)  and  signal 
acquisition  are  major  issues  in  the  design  of  LF  communications 
systems.  Consequently,  we  also  considered  these  topics  briefly 
in  Chapter  4  and  developed  appropriate  sequence  design  criteria. 


Once  our  sequence  design  strategy  was  clear,  we  proceeded  to 
consider  actual  designs  based  on  additive  and  multiplicative 
characters  in  Sections  4.6  and  4.7  respectively.  In  both 
sections,  we  began  by  reviewing  appropriate  definitions  and 
theorems.  Most  importantly,  we  reviewed  theorems  which  allowed 
us  to  bound  full  sums  of  additive  and  multiplicative  characters. 
These  full  sum  theorems  along  with  the  partial  sum  theorem 
discussed  above  allowed  us  to  bound  the  multiple  access,  skywave 
and  signal  acquisition  performance  of  our  sequence  sets.  In  both 
sections,  we  also  presented  a  brief  history  of  the  use  of 
additive  and  multiplicative  characters  for  code  or  sequence  set 
design.  Finally,  we  presented  an  additive  and  multiplicative 
character  sequence  set  design  for  our  LF/CDMA  communication 
system.  The  main  difference  between  the  two  sequence  designs  may 
be  the  number  of  phases  used  to  modulate  the  user  transmissions. 
The  additive  character  design  accommodates  N  users,  but  requires 
N  phases,  where  N  is  the  sequence  length.  The  multiplicative 
character  designs  can  use  many  fewer  than  N  phases,  but 
consequently  accommodates  many  fewer  than  N  users. 


APPENDIX  A 


DERIVATION  OF  THE  VARIANCE  OF 
MULTIPLE  ACCESS  INTERFERENCE 


A. I  Introduction 

In  this  appendix,  we  will  derive  the  variance  (or 
equivalently  the  second  moment)  of  the  multiple  access 
interference  found  in  our  r-phase/orthogonal  signal  set  CDMA 
system.  We  have  already  stated  the  result  in  Equation  (3.37)  and 
used  it  to  derive  many  important  insights  in  Chapter  3. 
Mathematically,  we  seek  an  expression  for 
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and  our  search  is  broken  into  the  following  sections. 

•  Define  convenient  symbology. 

•  Show  relationship  between  variance  of  real  interference 
and  complex  interference. 

•  Find  expression  for  variance  of  complex  interference. 

•  Find  expression  for  variance  of  real  interference. 


A. 2  Definitions  and  Symbols 

An  expression  for  l^a'1  is  given  by  Equation  (3.25)  and  that 
expression  is  repeated  below  except  that  the  notation  is 
streamlined  for  simplicity 

lx  =  &(c  Re{HyU)}  -  s  Im(HyU)})  (A-2) 


+  JsR(c  Re  {Hy  (  £+1)  }  -  s  Im{Hy  ( .1+1)  } ) 


where  x^=y  is  the  /Ith  sequence  and  x^=y  is  the  Ath  sequence  and 
where 


i4r(T(6)) 

R»Rr(T(S)) 

( 8 ) 

C=COSO)  T  vo/ 

c 

s=sina) 


c 

V*)=HxA 


u,v 


(A) 


(A. 3) 


An  expression  for  complex  multiple  access  interference  is  now 
defined  which  uses  similarly  streamlined  notation. 

e  H,  U)  +  R  e  H  ( 2.+1) 


(A. 4) 


where 


e=exp j^cT 


(6) 


(A. 5) 


A. 3  Relationships  Between  Real  and  Complex  Interference 

By  manipulating  Expressions  (A. 2)  and  (A. 4),  the  following 
relationships  can  be  derived 

Ix=JsRe{Ix}  (A. 6) 

|IXI  £  H\lx\  (A. 7) 

Eb{  I  XA  *  2  ^  £^tEb{lIA|2j  for  all  t  (A. 8) 

Additionally,  for  systems  where  the  carrier  frequency  varies 
quickly  when  compared  to  the  chip  waveform 

VIV1*!2’  "  8et"W2>  ia-9> 


The  last  equality  relates  the  second  moment  of  complex  multiple 
access  interference  to  the  second  moment  of  real  interference. 
It  is  the  most  important  of  the  above  expressions  and  will  be 
further  developed  in  the  next  section. 

A. 4  Second  Moment  of  Complex  Interference 

In  this  section  of  the  appendix,  we  seek  an  expression  for 

j 


w 


) 


I 


the  second  moment  of  the  complex  multiple  access  interference 


Vx'iv1*1  1 


This  search  is  made  convenient  by  proving  the  following  two 
lemmas . 


First  of  all,  we  wish  to  show  that 


(VM)  u^C{3)  UkUi  ~  6k-i 


(A. 10) 


where  uk  and  Uj  are  elements  of  the  vector  u,  which  belongs  to 
an  orthogonal  code  C$K  Consider  that  an  orthogonal  code  is  a 
set  of  vectors  where 


(1/N) 


N— 1 

z 

n=0 


* 

u  v  = 
n  n 


0  u/v 

1  u=v 


(A. 11) 


where  un  and  vn  are  vector  elements.  Equation  (A. 11)  may  be 
written  in  matrix  form  as  follows 

BB*T=NI  (A. 12) 


where  B  is  a  matrix  whose  rows  are  the  codewords  and  I  is  the 
identity  matrix.  The  left  hand  side  of  Equation  (A. 10)  may  be 
rewritten 


m 

=  B 


* 


(1/M)  T 


/  n  \ 


* 

U,  U. 


(A. 13) 


but  clearly 


(1/M)  I  uvut  =  BXB  =  (BB*T)T  =  NI  =  MI 

ueC  P  K  1 


which  proves  our  assertion. 

The  second  lemma  we  wish  to  prove  is  that 


(A. 14) 


E,  {H  U+a'  )H  (S,+a)  }  - 
b  x  y 


C  (a-a')  a  <  a' 
y,x  '  — 


CX,y {a'~a)  a,>  a 


(A. 15) 


We  proceed  as  follows 


Eb{HxU+a.')H*U+a)  }« 


(A. 16) 


Eb{  (Cx,uU+a,)+Cv,x(N"il'a,))  (Cy,uU+a)+Cv,y(N“*“a)  >  } 

“(l/M)  X  c;ullt,.|C  [««]  +  (l/M2)  X  I  C*  ( t+a '  )  C*  It,- 1- a) 

u  v  '  v,y 

+(1/M2)  lie  (N-4-a')Cv  U+a) 

u  v  ,x 

+  (l/M)  I  CV;X<N-,(-a')c;_ym-.«-a) 


-(l/M)  I  Xk-£-b  ^  y,  0  -  I  u,u*  +  0  +  0 

k=Jl+b  k  *  b  i=£+a  1  *  a  u  k  x 


4+b-l 

+ ( l/M)  I  x 


k+N-Z-b 


A+a-1 

I  y. 


i+N-2-a 


I  v.  v. 


If  we  substitute  in  our  result  from  the  previous  lemma,  which  is 
stated  in  Equation  (A. 10),  we  may  derive 


(l+a')Hy>U(V(I+a)  } 

(A. 17) 

N-l 

1+b-l  * 

ii*+a  yi-*-aXi-^-b  + 

xi+N-l-byi+N-l-a 

a  >  b 

N-l  * 

I  xi-£-hyi-l-a  + 

i=l+b  1  *  b  1  *  a 

1+a-l  . 

/Z0  yi+N-l-axi+N-l-b 

b  >  a 

cy,x(a-a,) 

a  ^  a' 

Cx,y(a,~a) 

a '  >  a 

This  final  equation  is  the  sought  after  result  for  the  second 
lemma. 

With  the  aid  of  Equation  (A. 15),  we  may  now  conveniently 
derive  an  expression  for  the  second  moment  of  the  complex 
interference.  We  begin  as  follows 

EbETlIV1X|2>  <A.18> 

N-l 

-e.  (i/T)  z  (2fr(|H  u)  r+jH  u+i)  r+|H  (z)  r+|H  u+i)  r- 
D  £=o  1  x  x  y  y 

2Re(Hx(Z)H*(SL)+Hx(Z+l)H*(Z+l)  }) 

f  (Hx  (  i+1)  Hx  U)  +Hx  (  i+1)  Hx  ( Z)  +Hy  ( " +1)  Hy  ( l) 

H*  (  1  +  1)  H  (  1)  +  H  ( l)  H*  ( 1+1)  +  H  ( 1  +  1)  H*  (  l)  )  ) 

/  y  x  y  x  y 


Applying  the  result  of  our  second  lemma  (Equation  (A. 15)),  we 
obtain 


EbET<IV1A|2> 

N-l 

=  ( 1/T)  I  0^T  4N  + 

1=0 

(2Retcx,x(1)  >+2Be{cy,y(i>  }<,y(-iH-c*,y(D )  ) 

=  (1/Tc)  r4N  + 

^M2RetCx,x(1)  >  +  2ReICy,y  <i;  V~1)+Cy(1>  ’  > 

A. 5  Second  Moment  of  Real  Multiple  Access  Interference 

By  applying  Equation  (A. 9)  to  Equation  (A. 19),  we  obtain  the 
following  result,  which  we  use  in  Chapter  3 

EbET((Iu"IX)2}  =  d/8Tc)(^r4N  +  (A. 20) 

?^p(2Re{C  (1)  }+2Re{C  „  (1)  }+C*  ..  (-1)+C*  (1))) 
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APPENDIX  B 


DERIVATION  OF  CHARACTERI STIC  FUNCTION  OF 


MULTIPLE  ACCESS  INTERFERENCE 


The  character  i  Stic  function  of  is  sought  where  l^'1 
and  1^0)  are  the  interference  terms  introduced  into  the  outputs 
of  the  and  correlators  by  the  (3^  interfering  user.  We 
denote 


I  <B)-I  -I,  = 

y  A  pi  A 

5ft(T)  Re(expju)  t(C*  _x  U(^)+Cv  x  (N-*))}  + 

y  A  y  a 

JsR(t)  Reiexpjoi  t  (C*  (£+l)+C  _  (N-S.-1)  )  } 

c  xu  x,,u  v,xy  xx 


4>(S)  (co)=E^B)E^6)E(i6)E(iCt)  {expjW(Iu-Ix)  } 


(B.l) 


(B .  2 ) 


The  above  expressions  may  be  combined  and  manipulated  to  provide 
the  following  key  relation 


„(6>(U|=E<e»E<8>E'8»E'“»( 
T  b  c  c 

N-l 

7 T  expjw  h^(x)  Re{|& 
n=l  1 


(B .  3 ) 


*  * 

—  O  a 

y,n-&  A  ,  n-2- 1  •  n-X. ' 


1  i  *  *  i 

7T  expjoi  ^(r)  Re{ \k  ...  „ ...  ,  a  ...  .  } 
_  1  u,n+N-£  A,n+N-£'  n+N- 1 

n=0 


N-l  , 
7T  expju)  ^R(t)  Re{|x 
n=H+l 


y ,  n-  9.-1  XA ,  n- 1-1  ^  an- i-1 ' 


7 T  expjoj  - 2 R  ( t  )  Re{|x 
n~6 


y  , n+N- ?.-l  XA, n+N- 1-1  an+N-l-l’ 


where  the  an  are  products  as  follows 


*  (a.)  ( S ) 

a  0=expjco  t  expjS  .  c  c  u 

n-x.  c  n-l  n-i 


n  n 


( B .  4  ) 


an+N-rexPi“oT  ^ 


(B.5) 


.  .  ~  *  (  V  \  (  Q  )  A 

>n-J-rexP3"oT  ^Vn-l  cnVl  un 


(B.6) 


(6) 


an+N-t-l=expj“cT  exp30n+N-l-l  cn+N-i-l  cn  vn  (B‘7) 


The  and  cj^  are  the  random  elements  of  the  coset  leaders 


for  the  ath  and  /3th  users  respectively,  hence 


XU  ,  n  Xy,nCn 


(a) 


(B .  8 ) 


u(8)c(8) 
N-l  N-l 


;(B)C(3) 
0  0 


(3)C(S), 

N-l  N-l; 


v 


(B .  9 ) 


The  $n  are  the  phase  angles  of  the  differences  Xju.,n_*X,n*  The  an 
are  complex  random  variables  with  unity  magnitude  and  phase  which 


112 


0 


is  uniformly  distributed  on  [0,2”).  Most  importantly,  they  are 
independent  for  all  n,  hence  we  may  derive 


M  1  ^  /\  /\  /\ 

b(3)(w)=E„  77  ( 1/2TT )  f  expjw^R(-)  |x  -x,  j 

1  n=0  ^  U,n  A,n 


jR  (- )  x,  -x,  |  cos8d5 


(B. 10) 


*  TT  ( 1/ 2 tt )J  expju)-iR(i)  |  x^ n-x^ n  |  cos0d8 
n=0  Q 


-(l/Tc)CC  V  J0(>i^(t)|lUin-lx,nl) 

"o  n  u 


w  A  > 

U  ,  n  A  ,  n 


TT  J0(^R(T)  |X  -XX  |)dT 
n— o 


This  expression  may  be  further  developed  by  considering 


Xn  n-x*  n | =/T/l-COS (0  -e,  _) 

U  /  n  a  ,  n  u ,  n  A  ^  n 


(B.ll) 


For  r-phase  orthogonal  codes,  codewords  are  orthogonal,  because 
the  set  of  phase  angle  differences  (0^  n”^X,  n^  are  evenly  spaced 
over  the  interval  [0,2JT).  Consequently,  we  may  write 


(  Q)  /■—  C  1  _ 

1  (w)-(l/T  )  TT  J  (»swR(t)  /2/1-cos  (2Trn/fl)  ) 

C  J  n  o 

0  n=0 
N-l 

o  77  J  ( '2 R  (x  )  /T/l-cos  (  2Tn/N)  )  dt 
n=0  ° 


( B . 12) 


This  is  the  expression  we  use  in  Section  3.4.3. 


APPENDIX  C 


CHARACTERISTIC  FUNCTION  APPROACH  FOR  DETERMINISTIC  CODES 

In  the  subsections  3.4.3  and  3.4.5,  we  have  used 
characteristic  functions  to  evaluate  the  probability  of  error  in 
systems  using  some  variety  of  random  coding.  This  technique  has 
resulted  in  probability  expressions  that  are  very  easy  to 
evaluate  numerically.  However,  Geraniotis  and  Pursley  (1982) 
have  also  used  the  characteristic  function  approach  with 
deterministic  sequences  in  their  antipodal  signalling  model.  In 
fact,  they  have  produced  algorithms  which  calculate  the 
probability  of  error,  and  the  computational  requirements  of  these 
algorithms  only  grow  linearly  with  K  or  N.  Unfortunately, 
similar  algorithms  for  our  model  are  much  less  efficient.  To  see 
why,  consider  the  multiple  access  characteristic  function 

<t>(SW=E.|6)E£6)  {exp (- jw  ( 2/N£p)  (I^S)  -l{B)  )  )  }  (C>1) 

As  shown,  this  function  depends  on  | u  and  X  and  must  be  averaged 
over  the  possible  sequence  overlaps  b^K  Consequently,  our 
algorithm  has  a  computational  requirement  proportional  to  KN^. 
Even  though  this  may  be  reduced  to  KN4  by  judicious  program 
design,  the  run  time  of  this  algorithm  is  prohibitive  for  all  but 
the  smallest  sequence  lengths. 

We  did  not  use  our  algorithm  to  investigate  average  system 
performance,  but  we  did  use  it  to  provide  some  results  on  the 
worst  case  probability  of  error  as  a  function  of  user  delays. 


This  was  done  for  the  following  three  reasons,  which  we  will 
expand  on  below. 


«  Worst  case  delay  performance  is  of  particular  interest 
for  the  LF  channel. 

•  The  random  coding  results  given  earlier  will  provide 
good  estimates  of  the  average  performance  of 
deterministic  sets.  However,  they  are  not  useful  for 
estimating  the  worst  case  delay  performance  of 
deterministic  sets. 

•  For  worst  case  performance,  our  algorithm  can  be 
made  computationally  more  efficient  after  making 
some  reasonable  assumptions. 

Worst  case  delay  performance  is  of  interest,  because 
unfavorable  relative  signal  delays  may  exist  for  some  LF  band 
users  indefinitely.  The  relative  signal  delay  between  two 
signals  will  change  if  the  delay  due  to  signal  propagation  or 
transmitter  timing  is  a  function  of  time.  However,  LF 
communications  usually  employ  a  groundwave  propagation  mode  which 
introduces  negligible  delay  change  as  a  function  of  time. 
Additionally,  the  transmitter  timing  must  be  rather  well 
controlled  to  allow  receiver  synchronization,  and  this  control 
prevents  rapid  relative  signal  delay  change.  This  latter  fact  is 
especially  true  if  the  LF  signals  are  also  serving  a  navigation 
role,  which  implies  extremely  stable  control  of  transmitter 
timing. 

The  random  coding  approach  employed  at  the  end  of  Chapter  two 


is  not  a  good  one  for  investigating  the  worst  case  probabilities 
as  a  function  of  user  delays.  For  the  random  coding  approach, 
the  maximum  probability  of  error  occurs  any  time  the  user  delays 
are  equal  to  an  integer  times  the  chip  width.  For  actual  codes, 
the  probability  of  error  for  any  set  of  user  delays  is  a  function 
of  the  various  partial  cross  correlation  functions. 
Consequently,  if  the  user  delays  were  constrained  to  integer 
multiples  of  the  chip  width,  the  resulting  probability  of  error 
would  still  vary  considerably.  However,  the  random  coding 
approach  does  provide  good  estimates  when  the  probability  of 
error  is  averaged  with  respect  to  the  user  delays. 

For  worst  case  performance,  our  algorithm  can  be  made 
somewhat  more  efficient,  provided  we  are  willing  to  make  some 
reasonable  assumptions: 

•  The  greatest  probability  of  error  occurs  when  all 

user  delays  are  integer  multiples  of  a  chip  width. 

•  The  worst  case  T ^ )  is  independent  of  tW)  for  all  (3 

net  equal  to  7. 

•  The  worst  case  delays  are  independent  of  eb/N0. 

An  algorithm  employing  these  assumptions  was  designed  and 

used  to  generate  the  upper  curves  in  Figures  C.l  and  C.2.  These 

curves  are  the  worst  case  probability  of  error  as  a  function  of 

user  delays  (max  Pr(e,T))  for  the  length  17  additive  character 
T 

codes.  The  lower  curves  are  the  average  probability  of  errors 
given  by  Equations  (3.58)  and  (3.59).  As  shown  in  the  figures, 
the  worst  case  probability  of  error  can  be  orders  of  magnitude 
greater  than  the  average  probability  of  error.  However,  Pr(e,~) 
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is  near  the  worst  case  curve  only  for  a  very  small  range  of  user 
delays.  Additionally,  a  small  set  of  results  indicates  that  the 
average  performance  of  the  additive  character  sets  is  very  close 
to  the  lower  curves.  Consequently,  the  average  performance 
curves  from  Chapter  3  are  good  design  tools  whenever  the  user 
delays  vary  even  slowly. 
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